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ABSTRACT 

New Chandra X-ray data and extensive optical spectroscopy, obtained with AAOmega on the 3.9 m 
Anglo-AustraUan Telescope, are used to study the complex merger taking place in the galaxy cluster 
Abell 2744. Combining our spectra with data from the literature provides a catalog of 1237 redshifts 
for extragalactic objects lying within 15' of the cluster center. From these, we confirm 343 cluster 
members projected within 3 Mpc of the cluster center. Combining positions and velocities, we identify 
two major substructures, corresponding to the remnants of two major subclusters. The new data are 
consistent with a post core passage, major merger taking place along an axis that is tilted well out of 
the plane of the sky, together with an interloping minor merger. Supporting this interpretation, the 
new X-ray data reveal enriched, low entropy gas from the core of the approaching, major subcluster, 
lying ~ 2' north of the cluster center, and a shock front to the southeast of the previously known 
bright, compact core associated with the receding subcluster. The X-ray morphology of the compact 
core is consistent with a Bullet-like cluster viewed from within ~ 45° of the merger axis. An X-ray 
peak ^ 3' northwest of the cluster center, with an associated cold front to the northeast and a trail 
of low entropy gas to the south, is interpreted as the remnant of an interloping minor merger taking 
place roughly in the plane of the sky. We infer approximate paths for the three merging components. 

Subject headings: galaxies: clusters: individual (Abell 2744) — X-rays: galaxies: clusters 



1. INTRODUCTION 

In a universe where structure grows hierarchically, clus- 
ters of galaxies are the latest structures to collapse and 
virialize (e.g., Springel et al. 2006). There are three main 
modes of mass accretion onto rich clusters of galaxies: 
the steady infall of matter from the surrounding fila- 
mentary large scale structures, the discrete accretion of 
group-sized objects, and the extreme event of a major 
cluster-cluster merger. The latter are the most ener- 
getic events known in the Universe (Markevitch et al. 
1998) and result in the violent reassembly of the cluster. 
Such a dramatic reconfiguration of the cluster results in 
a rapid change in the environment of its member galax- 
ies, although what effect this has the galaxies themselves 
is yet to be fully understood. This is in part due to the 
complex nature of cluster mergers and, hence, the diffi- 
culty in obtaining a detailed picture of their properties. 

The well known Butcher-Oemler effect (hereafter BO- 
effect; Butcher & Oemler 1978, 1984), that the fraction 
of blue galaxies in the cores of rich clusters is significantly 
lower at the present epoch than it was ~2.5 Gyrs or more 
ago, reveals that a significant fraction of cluster galax- 
ies have undergone rapid transformation in star-forming 
properties over the intervening period. Furthermore, the 
observed increase with redshift in the fraction of spiral 
galaxies in clusters at the expense of a commensurate 
decline in the SO fraction (Dressier et al. 1997; Fasano 
et al. 2000; Desai et al. 2007; Just et al. 2010) reveals 
a corresponding rapid evolution in galaxy morphology. 
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While the BO effect appears to be widespread, the scat- 
ter in the blue galaxy fraction is large at all redshifts 
(z > 0.2) and exceeds the uncertainties in the measure- 
ments (Butcher & Oemler 1984). This indicates some 
internal cluster-specific mechanism is responsible for the 
scatter and Margoniner et al. (2001) found that a large 
portion of this scatter can be attributed to the richness 
of a cluster, where less rich clusters have a higher blue 
fraction. Another intriguing possibility is that the scat- 
ter is driven by the hierarchical formation of clusters, 
in particular cluster mergers, which become more com- 
mon with increasing redshift. Here it is quite conceivable 
that they would lead to an increased blue galaxy fraction 
through triggering star formation in the member galax- 
ies (Kauffmann 1995; Metevier et al. 2000; Miller et al. 
2006). 

In this context, spectroscopic observations of the Coma 
cluster have revealed a tantalizing correlation between 
the spatial distribution of post-starburst galaxies and re- 
gions involved in merger activity (Caldwell et al. 1993; 
Caldwell & Rose 1997; Poggianti et al. 2004). These ob- 
servations strongly suggest that the merging process can 
affect the star-forming properties of the cluster galaxies. 
This spurred Caldwell & Rose (1997) to conduct a further 
study to search for "abnormal" spectrum galaxies, simi- 
lar to those BO galaxies found at higher redshift, in four 
additional nearby clusters, three of which were selected 
on the basis of harboring clear substructure. Although 
their study was limited to galaxies of early-type mor- 
phology, the results indicated that there is a significant 
fraction of abnormal-type galaxies in nearby rich clusters 
and, most intriguingly, the triggering of starbursts lead- 
ing to abnormal spectra occurs during the core passage 
phase of a merger. Recently, Hwang & Lee (2009) have 
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provided further evidence for this scenario by comparing 
the galaxy properties in the post-eore passage merger 
Abell 168 and the pre-core passage merger Abell 1750. 
Furthermore, radio observations have revealed that clus- 
ters which harbor evidence for major merger activity 
show an increase in radio activity amongst their galaxies 
(Miller & Owen 2003; Miller et al. 2006; Venturi et al. 
2001, 2000; Johnston-HoUitt et al. 2008). However, there 
are counter-examples (Venturi et al. 2002), suggesting 
that the merger phase is an important factor. According 
to the simulations of Bekki et al. (2010), this evolution 
may be driven by the significant increase in ICM pres- 
sure that a galaxy is exposed to during a cluster merger, 
in particular when the galaxy passes through regions af- 
fected by shocks during the core passage stage of the 
merger (see also Roettiger et al. 1996). Therefore, an 
understanding of the merger dynamics and merging his- 
tory of a cluster is critical when attempting to disentan- 
gle the effects leading to the triggering and halting of star 
formation in cluster galaxies. 

The rich X-ray luminous cluster Abell 2744 at z = 0.3 
(also known as AC118; Couch & Newell 1984) is an inter- 
esting test case, being both a well known merging cluster 
and one that exhibits a significant BO effect. Two pre- 
vious studies have addressed the dynamical state of the 
merger occurring in the core of Abell 2744. The first was 
the C/iandra-based study of Kempner & David (2004, 
hereafter KD04), who showed that Abell 2744 is under- 
going a major merger in approximately the north-south 
direction. The second was conducted by Boschin et al. 
(2006) who used 85 cluster member spectra to show that 
the merger has a significant velocity component along 
the line of sight, as suggested by KD04 (see also; Gi- 
rardi & Mezzetti 2001), and estimated the mass ratio 
of the merging subclusters to be ~ 3 : 1. This merger 
also quite naturally explains why Abell 2744 hosts one 
of the most luminous known radio haloes covering its 
central 1.8 Mpc, as well as a large radio relic at a dis- 
tance of about 2 Mpc from the cluster center (Giovan- 
nini et al. 1999; Govoni et al. 2001a,b). The presence 
of a significant blue galaxy excess — a blue fraction that 
was 2.2 ± 0.3 times that seen in the same core regions of 
nearby clusters was first measured photometrically and 
confirmed spectroscopically by Couch & Sharpies (1987), 
who found the cluster to have a blue galaxy fraction of 
~ 25%. Moreover, the spectroscopic study showed that 
the blue galaxy population in Abell 2744 was dominated 
by starburst and post-starburst galaxies. Abell 2744 
therefore provides an excellent laboratory for studying 
the link between major merger activity and star forma- 
tion activity in cluster galaxies. This paper reports the 
first step of this study, which is to determine more pre- 
cisely the phase and history of Abell 2744's merger, in 
particular whether it is in a pre-core or post-core pas- 
sage phase. 

As has been shown previously (e.g., Zabludoff & Zarit- 
sky 1995; Barrena et al. 2007; Maurogordato et al. 2008; 
Owers et al. 2009c; Ma et al. 2009, 2010), the combi- 
nation of high quality X-ray imaging and spectroscopy 
with comprehensive optical spectroscopy provides a pow- 
erful toolkit for disentangling the complex dynamics of 
cluster mergers. In this paper, we combine new Chan- 
dra observations, which provide an additional lOOks 
of data to the existing 25ks analyzed in KD04, with 



new AAT/AAOmega optical multi-object spectroscopy 
(MOS), which roughly doubles the number of spectro- 
scopically confirmed cluster members over a 3 Mpc ra- 
dius region, and also doubles the number of cluster mem- 
ber spectra within the central regions studied by Boschin 
et al. (2006). These data sets are used both to detect sub- 
structure in Abell 2744, which is related to the merger 
activity, and to provide a more up-to-date interpretation 
of the merging history. 

The outline of this paper is as follows. In Section 2 we 
present the reduction and analysis of the new Chandra 
data. In Section 3 present details of the MOS observa- 
tions and data reduction, the rcdshift catalog, and the 
cluster mcanbership allocation procedure. In Section 4 
we present the methods used for substructure detection. 
In Section 5 we present our interpretations of the nature 
of the structures detected in the optical and X-ray data. 
In Section 6 we present our merger scenario based on the 
observations and interpretations presented in the paper. 
Finally, in Section 7 we summarize our results. 

Throughout the paper, we assume a standard ACDM 
cosmology where Hq = 70kms~^, Qm = 0.3 and 
1 1.7. For the assumed cosmology and at the cluster red- 
shift (z=0.3064), l"=4.52kpc. 

2. CHANDRA DATA 

2.1. Observations and Data Reduction 

Abell 2744 has been observed with Chandra using both 
the ACIS-I and ACIS-S chip arrays and we summarize 
the dates, exposure times and Chandra ObsID's in Ta- 
ble 1. The data were reprocessed starting with the level 
1 event files and using CI AO version 3.4 with the latest 
gain and calibration files applied (CalDB version 3.5.5). 
All observations were performed in Very Faint (VF) data 
mode, allowing 5x5 pixel islands to be used in identi- 
fying cosmic ray events, significantly helping to reduce 
the particle background. Observation specific bad pixel 
files were produced and applied and events flagged with 
ASCA grades 1, 5 and 7 were excluded. The data were 
filtered for periods of anomalously high background as- 
sociated with flares. Standard procedures arc followed in 
doing this, i.e., for the ACIS-I observations, point sources 
and the region containing the diffuse cluster emission 
were excluded, and a light curve was extracted in the 
0.3 — lOkeV range binned in 259s intervals. The Sherpa 
Ic-clean.sl routine was then used to identify for removal 
time intervals when the count rate deviated by more than 
20% from the mean. No significant fiares were detected, 
and the cleaned exposure times are listed in Table 1. 

For background subtraction during spectral and imag- 
ing analyses, we use the blank sky observations^ appro- 
priate for the chip and epoch of the observation of in- 
terest. The background files were processed in the same 
manner as the observations, using the same background 
filtering, bad pixel and gain files, and were reprojected 
to match the observations. The backgrounds were nor- 
malized so that the source and background counts in the 
10 — 12keV range matched. We checked for excess soft 
Galactic X-ray emission in the observations by extract- 
ing spectra from source-free regions and comparing to 
spectra extracted from the background files, finding no 
significant difference. 

^ See http://cxc.harvard.edu/contrib/maxim/ELcisbg 
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Table 1 

Summajy of the Chandra observations. 



Array Obsid Date Exposure Time (ks) Cleaned Exposure Time (ks) 



ACIS-S 2212 2001 September 3 25.14 23.79 

ACIS-I 7712 2007 September 10 8.17 8.07 

ACIS-I 7915 2006 November 8 18.86 18.62 

ACIS-I 8477 2007 June 10 46.5 45.91 

ACIS-I 8557 2007 June 14 28.17 27.75 



2.2. Im,age Analysis 

The left panel of Figure 1 shows a background- 
subtracted, 0.5 — 7keV, exposure-corrected image of 
Abell 2744 which has been binned using the WVT bin- 
ning algorithm of Diehl & Statler (2006), which is a 
generalization of Cappellari & Copin (2003) 's Voronoi 
binning algorithm. The bin size is determined by the 
constraint that the signal-to-noise ratio must be ^ 5. 
Prior to binning, point sources detected with wavde- 
tect were removed and the regions wcirc; filled via in- 
terpolation from a surrounding background region us- 
ing dmfilth. The exposure map was generated using 
standard CI AO procedures^ and accounts for the effects 
of vignetting, quantum efficiency (QE), QE nonunifor- 
mity, bad pixels, dithering, and effective area. The en- 
ergy dependence of the effective area is accounted for 
by computing a weighted instrument map where the 
weights are computed from an absorbed MEKAL spec- 
tral model with Galactic absorption, temperature, red- 
shift and abundance set to the average cluster values ob- 
tained in Section 2.3.1. The right panel of Figure 1 shows 
the background subtracted, exposure corrected image of 
Abell 2744 with point sources included, and with a light 
Gaussian smooth (FWHM=4") applied. Abell 2744 is 
clearly a disturbed system, with a structure ^ 750 kpc 
to the northwest of the X-ray surface brightness peak 
(labelled northwestern interloper). Within the; main X- 
ray structure, there are a further three substructures: 
one 190 kpc south (labelled southern compact core), 
one ~ 480 kpc to the north (labelled northern core) and 
another ~ 220 kpc to the north of the X-ray surface 
brightness peak. The northern core furthest from the 
X-ray peak has a tail, which includes the second struc- 
ture north of the X-ray peak, extending southward and 
curving westward towards the peak in the X-ray emis- 
sion. The southern compact core has a blunted cone 
morphology with a prominent head and a fan-like tail of 
emission extending back towards the peak in the X-ray 
emission. There is also an edge in the surface brightness 
--^ 250 kpc to the southeast of the southern substructure. 

The northwestern interloper appears asymmetric, with 
an edge just east of north and an extended tail toward 
the south. There also appears to be a bridge of emission 
connecting it to the main cluster emission. Based on 
a shorter Chandra observation, KD04 tentatively identi- 
fied a cold front and bow shock on the eastern side of 
the northwestern interloper and also a fan-like feature 
extending towards the west. They concluded the north- 
western interloper was moving from west to east on its 
way towards the cluster core. However, the deeper ob- 
servations presented in Figure 1 do not confirm the exis- 

* see http://cxc.harvard.edU/ciao3.4/threads/expmap_acis_multi/ 



fence of a bow shock and indicate the cold front on the 
eastern side is part of a larger northern edge, meaning 
the interpretation of KD04 may not be the correct one. 

To emphasize the structures revealed in Figure 1 
we present the residual significance map, which high- 
lights low surface brightness structures, and the unsharp 
masked image, which is useful for revealing sharp struc- 
tures such shocks and cold fronts, in Figure 2. To pro- 
duce the significance map, we first fit the Chandra im- 
age with an azimuthally symmetric /3-model with best 
fitting values of Tc = 47i ± 8 kpc, j3 = 0.96 ± 0.13, 
R.A.=00'^14'°18.4^ and decl.=-30°23'14".0. During the 
fitting of the /3-model, the northwest substructure was 
excluded, along with point sources. This /3-model ac- 
counts for the main cluster emission and is subtracted 
from Chandra image to produce a residual image which is 
then smoothed with a Gaussian kernel with cr = 3". The 
resulting smoothed residual map is then divided by an 
error map to produce the residual significance map (see 
Neumann & Bohringer 1997; Owers et al. 2009c, for more 
details). The left panel in Figure 2 shows that the north- 
west interloper and its extended tail to the south are 
detected with high significance, while the northern core 
and southern compact core, along with the tail near the 
northern core, are also detected. The unsharp masked 
image is produced by dividing a background subtracted, 
exposure corrected image that was smoothed with a 4.4" 
Gaussian by the same image smoothed by 22", corre- 
sponding to physical scales of 20 kpc and 100 kpc at the 
cluster redshift. The unsharp masked image is shown in 
the right panel of Figure 2, where it can be seen that the 
structure in the central cluster region is very similar to 
that seen in the residual significance map. However, the 
unsharp masked image accentuates the edge to the north 
of the northwestern interloper and also a finger of emis- 
sion extending from the northwestern interloper towards 
the south. 

2.3. Thermodynamic maps 
2.3.1. Temperature maps 

To better determine the nature of the features seen in 

Figures 1 and 2 and to identify other regions affected 
by the merger activity, we produce two maps of the 
projected temperature distribution for Abell 2744 using 
complementary techniques. The first temperature map 
is produced using a similar method to that developed 
by O'SuUivan et al. (2005) and Maughan et al. (2006), 
which is outlined in full in Randall et al. (2008). Briefly, 
the X-ray image is binned to a 1".968 x 1".968 pixel size 
and at the position of each pixel, source and background 
spectra are extracted in a circular region surrounding the 
pixel. The radius of this region is defined adaptively so 
that the extracted region contains 1500 background sub- 
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Figure 1. Background-subtracted, exposure-corrected Chandra images. The image in the left panel is binned to S/N=5 using the WVT 
algorithm and the point sources have been removed and filled using dmfilth. The image in the right panel shows the image with point 
sources, binned to 1".968 pixels and smoothed with a Gaussian kernel with FWHM=4". 




30.0 25.0 20.0 15.0 10.0 05.0 0:14:00.0 
R.A. (J2000) 



30.0 25.0 20.0 15.0 10.0 05.0 0:14:00.0 
R.A. (J2000) 



Figure 2. Left panel: Residual significance image produced using a similar method as in Owers et al. (2009c), i.e., an azimuthally 
symmetric /3-model is fitted to the data (excluding the northwestern interloper and point sources) and then subtracted to produce a 
residual image. The residual image is then smoothed and divided by an error map to produce the residual significance map. Contours are 
spaced linearly in steps of la in the range 1 — ITa. Right parte/.Unsharp masked image created by dividing an exposure and background 
corrected image smoothed using a Gaussian kernel with a = 4.4" (20 kpc) by one smoothed by cr = 22" (100 kpc). 



tracted counts. The process is expedited by extracting 
the response files from a coarser grid. For each spec- 
trum, an absorbed MEKAL model is fitted in the energy 
range 0.6 — 9.5 keV, with the absorption column den- 
sity set at the Galactic value, Nu — 1.61 x 10^°cm~^ 
(Dickey & Lockman 1990), the abundance set to the 
global value of Zave — 0.262 (see below) and the tem- 
perature free to vary. This temperature map is pre- 
sented in the top left panel of Figure 3. We note that 
this method produces a temperature map with correlated 



pixels and the scale over which the pixels are correlated 
is not clear from the temperature map alone. There- 
fore, as a check we produced a second temperature map 
using the WVT algorithm to bin the image such that 
each (uncorrelated) spatial bin contains ^ 1500 back- 
ground subtracted 0.5— 7keV counts. We select a number 
of regions of interest from the temperature map shown 
in the top left panel of Figure 3 (marked with crosses) 
and we require the WVT binning algorithm to place a 
bin centered on each of these regions. For each bin we 
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extract observation-specific source spectra, correspond- 
ing count-wciglitcd responses, and background spectra 
from the blank sky observations. The background spec- 
tra had their exposure times normalized so that the 9.5- 
12KeV count rates match the observed rates. The spec- 
tra were re-binned so that there was at least 1 count 
per energy bin. The re-binned observation-specific spec- 
tra were fitted simultaneously over the 0.5 7 keV energy 
range within the XSPEC package (Arnaud 1996) using 
an absorbed MEKAL model (Kaastra 1992; Liedahl et al. 
1995) with abundance, column density and redshift fixed 
to those values used for the non-tessellated map. The 
Cash statistic was minimized during fitting. This tessel- 
lated temperature map is shown in the middle left panel 
of Figure 3. The tessellated and non-tessellated temper- 
ature maps are in excellent agreement. 

We have produced significance maps for both the tes- 
sellated and non-tessellated temperature maps which 
show how significant the spatial variations in temper- 
ature are when compared to the global average temper- 
ature, kTave- The significance, A(fcT), is defined at each 
pixel as 



A(fcT) 



(1) 

where fcTpix is the temperature map value for the pixel 
of interest, cr(fcTpix,hi) andCT(kTpix,io) are the upper and 
lower 68% confidence limits for fcTpix, respectively, and 
c7(A;Tave,hi) and cr(kTave,io) are the upper and lower 68% 
confidence limits for kTave, respectively. The significance 
maps are shown in the top right and middle right panels 
of Figure 3 for the non-tessellated and tessellated maps, 
respectively. The global average temperature, fcTavo, 
was determined from fitting spectra extracted from a 
~ 900 kpc radius region containing the majority of the 
cluster emission with an absorbed MEKAL model, as 
outlined above for the tessellated maps, with the excep- 
tion that the abundance was free to be fitted. We find 
kTa-vc = 9.07 ± 0.15 keV, while the global average abun- 
dance Zavc = 0.262 ± 0.025. These values are consistent 
with the values presented in Zhang et al. (2004). 

The temperature maps reveal a rich diversity of patchy 
temperature structures, the most significant of which in- 
clude: (i) A prominent cool ~ 5keV region associated 
with the northwest interloper where the coolest gas is 
coincident with the compressed X-ray isophotes, indicat- 
ing there is a cold front there. We also note a hint of a 
cool swirl heading to the south of the northwestern in- 
terloper and coincident with the asymmetry seen in the 
X-ray brightness (Figure 2). (ii) A cool ^ 7.5 keV region 
associated with the southern compact core which is sur- 
rounded by hot > 12 keV gas, the hottest of which resides 
just to the east of the substructure, (iii) A ~ 700 kpc 
ridge of hot > llkeV gas running towards the north 
which separates the northwestern interloper and north- 
ern core, (iv) While the northern core and its tail do 
not appear to have a significantly different temperature 
to the average temperature, they do appear to be signif- 
icantly cooler than the regions directly to the west and 
hotter than the region directly to the east. This is espe- 
cially clear when considering the tessellated map (middle 



left and right panels of Figure 3). (v) The peak in the 
X-ray surface brightness does not host the coolest gas, 
as would be expected in a relaxed, cool core cluster. 

2.3.2. Pseudo-Entropy, Pseudo-Pressure and Abundance 

maps 

From the non-tessellated temperature map, we have 
generated projected pseudo-pressure and pseudo-entropy 
maps, similar to those presented in Finoguenov et al. 

(2005). These are shown in the bottom left and bot- 
tom right panels of Figure 3, respectively. The pseudo- 

— 1/3 

entropy at each pixel is = fc^pix^pix > while the cor- 

1 /2 

responding pseudo-pressure is P = kXpixA^^^ where Api^ 
is the normalization determined during the fitting of the 
spectra corrected for exposure and normalized by the 
area of the extraction region. The northwestern inter- 
loper stands out as having the lowest pseudo-entropy in 
the cluster complex, while the swirl of cool gas, hinted at 
in the temperature maps, is clearly seen as low entropy 
gas extending towards the south. The non-hydrostatic 
nature of Abell 2744 is further highlighted by the asym- 
metric, patchy pressure distribution which peaks at the 
X-ray surface brightness peak and has two high pressure 
fingers trailing off to the southeast and northeast. 

Finally, we present an abundance map in the lower left 
panel of Figure 4. The technique outlined in Section 2.3.1 
used to generate the non-tessellated temperature map 
was again applied. However, the abundance was free 
to vary during the fitting while, in order to allow more 
robust abundance determinations, the spectra were ex- 
tracted from larger regions which contain ~ 3000 counts. 
The abundances were measured relative to the solar pho- 
tospheric values of Anders & Grevesse (1989). Included 
in the top left panel of Figure 4 is the corresponding 
temperature map which is in broad agreement with the 
temperature maps presented in the top and middle left 
panels of Figure 3. The most notable enhancements in 
abundance seen in the maps presented in the lower pan- 
els of Figure 4 are coincident with the southern compact 
core, the northern core and the northwestern interloper. 
We identify the northern core as the highest metallicity 
region within Abell 2744. Furthermore, there is an anti- 
correlation between the abundance and the temperature 
insomuch as the hottest regions also appear to have the 
lowest abundances. 

2.4. Regions of Interest 

We now explore the nature of several regions of inter- 
est which have been identified from Figures 1, 2 and 3. 

The selected regions are overlaid onto the temperature 
and abundance maps in Figure 5. For each region of in- 
terest the source spectra, along with the corresponding 
backgrounds and responses, are extracted and fitted in 
XSPEC with an absorbed MEKAL model. During the 
fitting, the column density is fixed to the Galactic value 
while the temperature, abundance, normalization, and, 
where the data are of a high enough quality, the redshift 
are allowed to free to vary during the fit. The best-fitting 
parameters for each of the regions are presented in Ta- 
ble 2. 

2.4.1. Northern Core and Surrounding Structures 
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Figure 3. Thermodynamics maps derived from the Chandra data. All maps have Chandra surface brightness contours overlaid. Top left: 
Temperature map using the method outlined in Randall et al. (2008). The cross points indicate regions where we require a WVT bin to be 
centered (see text). Middle left: Temperature map determined using the WVT binning algorithm to define regions. The color bars on the 
top and middle left panels show the temperature scale in keV. Top and middle right: Significance of the difference in temperature at each 
pixel compared to the global temperature where positive (negative) values indicate a higher (lower) temperature. The color bars show the 
significance scale which has units of the number of as deviation from the global temperature. Bottom left and right: Pseudo-entropy left 
and Pseudo-pressure right maps derived from the temperature and normalization. The units of the pseudo-entropy and -pressure maps are 
arbitrary. 
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Figure 4. Top left and right Same as the top left and right panels in Figure 3, but using an extraetion region which contains 3000 
background subtracted counts and allowing the abundance to vary during the fit. Bottom left: The abundance map where the color bar 
shows the corresponding abundance scale relative to solar. Bottom right: The significance of the difference of the abundance at each pixel 
relative to the global abundance measure. The color bar gives the scale in units of a. 
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Table 2 

Summary of fits to spectra extracted from the regions of interest. 



Region Description kT (IceV) Abundance (solar) Zgaa Net Counts 



Main Cluster 

Rl central QOOkpc 'd-^^tl'Xl 0.262l^;^25 O-SOSSlglooao 94206 



Northern Core Surrounding Structures 

NCI NC head 7.36155:50 0-480l°;J^° 0.3035i^;'*|^ 3103 

NC2 NC Tail 8.691°-^^ 0.424l°;J^^ 0.2846l!5:o?29 ^447 

NC3 structure south of NC 0.404l°;j3^ 0.2967lo g?i4 2962 



Southern Compact Core and Surrounding Structures 

SCCl SE cool core 7.70i;!;'^^ 0.392l°:J|^ 0.2950l°:|j'i'^^ 2059 

SCC2 Outsidc2 SE cool core 9.74l};^7 0.329to 2?3 0.3064 1551 

SCC3 Hot reg north of SE cool core 15.56li gJ O-OioiS 0.3064 2986 

SCC4 Inside edge 15.83t^;?Ej O-l^Stolag? 0.3064 1800 

SCC5 Outside edge 8.6ll J fj 0.607lo;2|? 0.3064 1260 



Northwestern Interloper and Surrounding Structures 

NWl NW inside front 6.33l°;fg 0.233l°;J^^ 0.308611;;°^^^ 1270 

NW2 NW outside front 10.45tj|Q OAOi^}]:'^'''^ 0.3064 1263 

NWl NW inside front deprojected 5.09lj];**^ 0.2931°;^^^ Q.ZO^bt!a^ouH 1270 

NW3 NWtail adjacent east 10.441J;™ 0.300t°;2oi 0.2935":S!.o^f43 1642 

NW4 NWtail 8-12lo:99 0.383l°;J^* 0.32150:{3io°/o7 1686 

NW5 NWtail adjacent West 8.28l};g^ OmQt°o\% 0.3064 1723 

NW4+NW5 Comb NW4 and NW5 8.34l°;^g 0.190l°;'^^";^*' 0.3189t°;°J^^ 1723+1686 

NW6 NWtail adjacent West 2 11.28t^ 37 0.007l°;^|^ 0.3064 1635 



Miscellaneous Regions 

MISCl ridge c 12.46l:i;i4 ^-^^"^-lAVi 0.3064 4258 

MISC2 SB peak 10.33t|3-|| O-SlSlo'ioi 0.3189l!5:;;'|'?f) 4793 

MISC3 Hot ridge north 13.14ti|^ 0-353lo:2oi 0.2813l!5:"i32 2470 

MISC4 low abundance 11.46lg;^^ 0.180lg:g^^ 0.3112l[j;™j'^^ 7897 



The abundance map (bottom left panel of Figure 4) 
reveals the northern core has significantly higher metal- 
licity than the surrounding gas, while the temperature 
difference is not significant. Furthermore, both the resid- 
ual significance and unsharp masked maps (Figure 2) re- 
veal a curved tail of emission trailing to the south of the 
substructure, as well as a second component ^ 270 kpc 
to the south. We extract and fit spectra from regions 
surrounding these structures (NCI, NC2 and NC3 in the 
right panel of Figure 5). We find that the northern sub- 
structure (NCI) does in fact have a significantly lower 
temperature than the global average, while the NC2 and 
NC3 regions have temperatures which are consistent with 
the global average within the quoted confidence levels. 
Also, we find that at the 68% confidence level all three re- 
gions harbor gas with higher abundances than the global 
value. 

2.4.2. Southern Compact Core and Surrounding Structures 

There are several regions of interest surrounding the 
southern compact core including the compact core itself 
(SCCl in the left panel of Figure 5) and two hot regions 
identified in the temperature maps which straddle the 
core on the northeastern and southwestern sides (SCC2 
and SCC3 in the left panel of Figure 5). While the de- 



tailed spectra extracted for the SCC2 region reveal the 
gas there is not significantly hotter than the global av- 
erage temperature, the region to the northeast (SCC3) 
does contain significantly hotter gas which also has a very 
low abundance. We confirm that the southern compact 
core harbors gas which is both cooler and of a higher 
abundance than the global averages. 

We attempt to quantify the strength of the surface 
brightness edge just southeast of the southern compact 
core (Figure 1) by fitting a spherically symmetric bro- 
ken power law density model centered at the center of 
curvature of the edge (using the same method as Owers 
ct al. 2009b). The resulting fit gives a jump in density 
at the edge of 1.601q;^^. The surface brightness profile 
and the best fitting broken powerlaw density model are 
shown in Figure 6. We note that in Obslds 7712 and 
7915, the chip gap lies very close to the position of the 
edge. Therefore we repeat the fit using only Obslds 8477 
and 8557. We find that the best fitting surface bright- 
ness model gives a density jump of 1.5llo:i55 consistent 
with the measurement obtained using the full dataset. 
The temperature measured just inside the edge (SCC4) is 

kT=15.83t2:f keV, while kT=8.6lL-f42 keV in the SCC5 
region just outside the edge. The change in pressure 
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Figure 5. Regions of interest labelled as in Table 2 shown overlaid on the temperature map {left panel) and abundance map {right panel). 

geometry assumed in the density model accurately repre- 
sents the physical geometry of the edge, the temperature 
just inside the edge will be artificially boosted by hot gas 
from outside the edge which lies along the line of sight. 
To account for this, we determine the deprojected tem- 
perature inside the edge by including a second MEKAL 
component in fitting the spectrum extracted from the 
NWIl region. This second MEKAL component has tem- 
perature, normalizations and abundance frozen to the 
values determined for the NWI2 region, with the normal- 
izations multiplied by a constant which accounts for the 
different emission measures due to the different volumes 
probed (see Owers et al. 2009b, for a detailed explanation 
of the procedure). The deprojected temperature in the 
NWIl region is 5.09^q'7^ and combining this with the 
measured density jump and the temperature measured 
in the NWI2 region, we find that the pressure across the 

front, given by ine,i„kTin)/ine,outkTout) 
continuous and, thus, the edge is a cold front 

We extract spectra for 4 regions running from cast to 
west across the northwestern interloper's low entropy tail 
(NWI3, NWI4, NWI5, NWI6 in the left panel of Fig- 
ure 5). The temperature map (top left panel in Figure 3) 
indicated that the low entropy tail may harbor cooler gas 
than the regions immediately to the east and west. Re- 
gion NWI4 coincides with the lowest projected entropy, 
while region NWI5 coincides with the lowest tempera- 
ture gas seen in the temperature maps in the vicinity of 
the low entropy tail. Regions NWI3 and NWI6 are the 
regions immediately east and west of the tail, respec- 
tively. The temperature and abundance measurements 
are listed in Table 2 and we find that the temperatures in 
NWI4 and NWI5 are indeed lower than those measured 
in NWI3 and NWI6, but only at the 68% confidence level. 



Figure 6. Surface brightness profile across a sector containing 
the edge to the southeast of the southern compact core (see re- 
gions SCC4 and SCC5 in Figure 5). The red curve shows the sur- 
face brightness profile for the best fitting broken powerlaw density 
model. 

across this edge is {ne,inkT^n)/{ne,outkTout) = 2.94l[];^|, 
which indicates a significant pressure jump and thus that 
the edge is likely a part of a shock front. 

2.4.3. Northwestern Interloper and Surrounding Structures 

As noted in Section 2.2, the northwestern interloper 
contains a sharp surface brightness edge to the north. 
The temperatures just inside and just outside the front 
(regions NWIl and NWI2 in the left panel of Figure 5) 
are 6.33^!^:?^ and lOAbtHl keV, respectively. Figure 7 
shows the surface brightness profile across the edge, 
which we have fitted with a spherically symmetric bro- 
ken power law density model centered at the center of 
curvature of the edge (using the same method as Owers 
et al. 2009b) with a best fitting density jump of2.13t°:?° 
(uncertainties are 68% confidence limits). Providing the 



i.04t°:i is 



2.4.4. Miscellaneous regions with high temperature and/or 
low abundance. 
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Figure 7. Surface brightness profile across a sector containing 
the NW substructure cold front (see regions NWIl and NWI2 in 
Figure 5. The red curve shows the surface brightness profile for 
the best fitting broken powerlaw density model. 

There are several regions in Figures 5 which have 
either high temperature (MISCl, MISC2 and MISC3) 
and and/or an anomalously low abundance (MISC4). 
The first of the hot regions (MISCl) is roughly coin- 
cident with a finger of excess emission clearly seen in the 
unsharp masked and residual significance images (Fig- 
ure 2) which extends ~ 450 kpc from the X-ray sur- 
face brightness peak and towards the north-northeast (la- 
belled Ridge C in KD04). The second region (MISC2) 
is coincident with the peak in the X-ray surface bright- 
ness while the third region (MISC3) is a streak of hot 
gas which extends ^ 700 kpc from the peak in the X-ray 
surface brightness towards the north. We also extract a 
spectrum for the region which encompasses the low abun- 
dance finger running from the peak in the X-ray surface 
brightness towards the north-northwest (MISC4). 

We find that the MISC2 region is marginally hotter 
than the cluster mean temperature and has an abun- 
dance which is statistically consistent with the mean 
value. The two other regions selected based on their high 
temperature, MISCl and MISC3 are both significantly 
hotter than the mean temperature, consistent with the 
temperature map results. The abundance of the MISC3 
region is consistent, within the uncertainties, with the 
mean cluster value, while the abundance of the MISCl 
region is low when compared to the cluster mean value. 
The MISC4 region, selected based on its low metallic- 
ity, does have a best fitting abundance which is lower 
than the mean cluster value, but this is not statistically 
significant. 

3. OPTICAL DATA 

In this section, we outline the new spectroscopic mea- 
surements taken using the AAOmega multiobject spec- 
trograph (MOS) which, when combined with the avail- 
able redshifts from the literature, provide a sample of 
343 spectroscopically confirmed cluster members which 
we utilize to detect substructures and to constrain the 
line-of-sight dynamics of the merger. 

3.1. Sample selection, Observations and Completeness 

Here we outline the selection of galaxies targeted for 
spectroscopic observations, the observation and reduc- 



tion of the spectra, the precision of our redshift measure- 
ments and the spectroscopic completeness of our sample. 

3.1.1. Photometric Catalog 

The primary photometric catalog used for selection of 
the spectroscopic targets was drawn from the Supercos- 
mos Sky Survey (SSS) server^ using the object catalog 
extraction function. We include only objects within a 
15' (- 4Mpc) radius centered at R.A.=00hl4™19^5, 
decl.=-30°23'19".0 (i.e. the center defined by Abell 
et al. 1989) and brighter than a magnitude of rF=21. 
We use the SSS classification fiag to remove those ob- 
jects with rp< 19 and not classified as galaxies (the SSS 
image classification is most reliable for rp < 19 Hambly 
et al. 2001) and do not apply any image-classification- 
based filtering for fainter objects. We cross-match our 
SSS catalog with the photometric catalog of Busarello 
et al. (2002) and include those objects which are not 
detected in the SSS, but appear in the Busarello et al. 
(2002) catalog and meet our magnitude limit. We ex- 
clude a number of objects which are not cluster members 
based on VLT VIMOS spectroscopy (Filiberto Braglia, 
priv. comm.; see also Braglia et al. 2009). The resulting 
catalog contains 1443 targets. 

3.1.2. AAOmega Observations 

Observations were taken during the nights of 2006 
September 12-18 using the AAOmega MOS on the 
Anglo-Australian Telescope (AAT; Saunders et al. 2004; 
Smith et al. 2004; Sharp et al. 2006). AAOmega consists 
of 392 2" diameter fibers available for target allocation 
across a two degree diameter field of view. During our 
observations, we used the medium resolution (R 1300) 
580V (blue arm) and 385R (red arm) gratings which pro- 
vide a wavelength coverage of 3700-8800A and spectral 
resolutions of 3.5A and 5.3A at 4800A and 7250A, re- 
spectively. 

The target density greatly exceeds the maximum fiber 
packing density, requiring multiple configurations in or- 
der to obtain a high level of spectroscopic complete- 
ness, particularly in the dense core regions of the clus- 
ter. We utilize the simulated annealing algorithm (Mis- 
zalski et al. 2006) option within the AAOmega spe- 
cific CONFIGURtP software to generate the fiber al- 
location files. This software operates most efficiently 
with an input catalog of ^ 500 — 800 objects, partic- 
ularly when configuring clustered fields. These require- 
ments are met by partitioning the ~ 1500 object pri- 
mary target catalog into subsamples based on the target 
magnitude. The magnitude bins used were rF< 19.63, 
19.63 <rF< 20.05, 20.05 <rF< 20.37, 20.37 <rF< 20.75 
and 20.75 <rF< 21.00. 

The limiting magnitudes of the subsamples serve two 
purposes: 1) they allow a fairly random subsampling 
of the densest regions, thereby reducing the computa- 
tion time for the CONFIGURE software and 2) allow 
a shorter observation time to reach the desired S/N 
limit for the brighter subsamples, enabling a more ef- 
ficient observing strategy. The limits are chosen to allow 
approximately two configurations per subsample, with- 
out under-sampling fiber allocations. The configurations 

^ See http://www-wfau.roe.ac.uk/sss/ 

^ See http: / /www. aao.gov.au/AAO/2df/aaomega/aaomcga_SAconfigure. html 
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are done in sequence, starting with the brightest mag- 
nitude subsample and continuing to fainter siibsamples. 
The unallocated objects in each subsample were propa- 
gated through to the next subsample in the configura- 
tion sequence. During the fiber configuration process, 
the CONFIGURE software attempts to allocate objects 
of higher priority ahead of lower priority objects, where 
prioritization is based on a user-defined rank assigned to 
each object. We required a high level of spectroscopic 
completeness in the central regions, so objects within 
500kpc of the cluster center were ranked highest. Rank- 
ings were decreased as a function of radius in 500 kpc 
radial bins out to a limit of 4 Mpc. 

This method produced nine configurations and we 
summarize the observations of these configurations in 
Table 3. For each configuration, ~ 40 fibers were al- 
located to blank sky regions for sky subtraction and a 
further 8 fibers were allocated to bright stars for acqui- 
sition and guiding. Differential atmospheric refraction 
limits the length of contiguous exposures for any one con- 
figuration to about three hours, therefore configurations 
requiring exposures longer than 3 hours were observed 
across several nights. For each configuration, we took a 
5s dome flat and 40s FcAr arc lamp exposures for the 
purpose of flat fielding and wavelength calibration of the 
spectra (Note: for configurations observed across several 
nights, night-specific calibrations were taken). The data 
were reduced using the AAO 2dFDR pipeline software'' 
which rcdiiccs the red and blue arm data separately, ex- 
tracting sky subtracted, flat fielded and wavelength cali- 
brated data for each exposure frame. The reduced frames 
are then co-added (using a weight derived from the to- 
tal fiux in each frame), at which point cosmic rays are 
identified and rejected. After co-addition, the red and 
blue-arm spectra are spliced together. At this point, the 
data wore processed using an IDL code which reduces the 
residuals left from sky subtraction via a Principal Com- 
ponent Analysis (PCA) process which is adapted from 
software developed by Wild & Hewett (2005, Emily Wis- 
nioski priv. comm. See Drinkwater et al. 2010 for a more 
detailed explanation). 

The reduced, combined, PCA corrected data are then 
assigned redshifts using the RUNZ code, which was writ- 
ten by Will Sutherland for the 2dF Galaxy Redshift Sur- 
vey (2dFGRS Colless et al. 2001). The RUNZ code has 
been adapted for use with AAOmega spectra and sig- 
nificantly enhanced from its earlier versions. In particu- 
lar, emission line redshifts measurements have been sig- 
nificantly improved for the WigglcZ siirvey (Drinkwater 
et al. 2010). The code measures redshifts using two meth- 
ods: 1) for absorption line spectra the cross-correlation 
method of Tonry & Davis (1979) is employed with a suit- 
able library of template spectra, and 2) for emission line 
galaxies the redshifts are measured for each line by fit- 
ting a Gaussian, with the final emission line redshift be- 
ing determined from the variance weighted mean of the 
redshifts determined for each emission line. The quality 
of the assigned redshift was determined by visual inspec- 
tion and an integer value classification, Q, was assigned, 
where Q ranges from 1 to 6, rising the same scheme 
adopted for the 2dFGRS, as described in Owers et al. 
(2009a, c). Briefiy, objects with Q < 2 are assigned no 

See http://www.aax3.gov.au/AAO/2df/aaomega/aax3mega_softv 



redshift or an unreliable redshift, objects with Q = 3, 4 or 

5 have reliable redshifts (with Q —5 having a template- 
like quality spectrum) and objects with Q =6 are stars 
or non-extragalactic objects. The AAOmega observa- 
tions produced reliable {Q =3,4 or 5) redshifts for 983 
extragalactic objects and 276 stars. 

3.1.3. Redshift uncertainties 

Of the 983 Q = 3, 4 or 5 AAOmega observations, there 

were 38 galaxy spectra with repeat observations. These 
repeat observations were used to check the uncertainties 
assigned by RUNZ to redshift measurements. We use a 
robust biweight estimator to determine the mean and dis- 
persion of the distribution of redshift differences, where 
the redshift measurement with the lower quality (defined 
by ether its lower Q value, or higher error measurement) 
is subtracted from the higher quality measurement. We 
find a mean difference of Acz = 1 ± 13 kms^^ and a dis- 
persion of 87 ± 21 kms~"^ from which wc can infer a sin- 
gle redshift measurement uncertainty of 61 ± 15kms~^. 
The median value of the individual redshift uncertainties 
assigned by RUNZ to non-stellar objects was 66kms~'^ 
(after excluding 6 outliers with error values greater than 
1000 km s~^), which is consistent with the uncertainty 
derived from the repeat observations. 

Independent checks on the redshifts and their uncer- 
tainties can be made by cross-correlating the AAOmega 
catalog with the Couch & Sharpies (1987), Couch et al. 
(1998), Boschin et al. (2006) and Braglia et al. (2009) 
redshift catalogs. Further to these catalogs, we add 
archival 2dF data, which is reduced using 2dFDR and 
redshifted using the RVSAO package in IRAF (Kurtz 

6 Mink 1998). The results of the comparisons are pre- 
sented in Table 4 where the columns are: Origin of exter- 
nal redshift catalog for comparison (Col. 1); number of 
redshifts in common with the AAOmega catalog, N match, 
where catalog members are deemed to be matches if their 
positions differ by less than 3" (Col. 2), the mean red- 
shift difference, Acz — czours — cz^xt., as determined by 
the biweight estimator (Col. 3), the dispersion in the 
Acz distribution, a{/S.cz), determined with the biweight 
estimator (Col. 4), the redshift uncertainty for the ex- 
ternal redshift measurements (Col. 5), the quadrature 
difference between (t(Acz) and the external redshift er- 
ror which gives an external estimate of the uncertainty 
on the AAOmega measurements (Col. 6) and the ratio 
of the external uncertainty measurement to the internal 
uncertainty measurement derived above (Col. 7). We 
note that the results in Table 4 show our internal red- 
shift uncertainties are systematically underestimated by 
a factor of ~ 1.8. Wc conclude that the redshift precision 
of the AAOmega measurements are ^ llOkms"^. 

3.1.4. Combined Redshift Catalog 

As can be seen from Table 4, there are significant off- 
sets from a mean of zero for all of the Acz distributions, 
aside from the comparison with the Boschin et al. results. 
In order to ensure the external redshift measurements 
are consistent with the AAOmega ones, prior to combin- 
ing the catalogs we add the pertinent value of Acz to 
each redshift measurement in the Braglia et al.. Couch 
et al. and Couch & Sharpies catalogs. When combin- 
'iH^lf^'feitalog, for consistency we always use AAOmega 
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Table 3 

Summajry of the AAOmega observations. 



Field 


Date 


rp Limit 


T 


^objects 


^Redshift 


Seeing 


1 


2006 Sep 12 


19.63 


5 X 1800s 


192 


180 


1".2 


2 


2006 Sep 13, 14 


19.63 


4 X 1800 + 3 X 1800s 


134 


126 


1".4 


3 


2006 Sep 12 


20.05 


10 X 1800s 


172 


161 


1.2-1".4 


4 


2006 Sep 13, 14, 15 


20.37 


4 X 1800 + 6 X 1800 + 5 x 1800s 


181 


147 


1.0-1".5 


5 


2006 Sep 13, 14, 15 


20.37 


1 X 1800 + 5 X 1800 + 4 X 1800s 


130 


114 


1.0-1". 5 


6 


2006 Sep 15, 16 


20.75 


5 X 1800 + 7 X 1800s 


190 


152 


1.0-1". 5 


7 


2006 Sep 16, 17 


20.75 


5 X 1800 + 5 X 1800s 


133 


109 


l."0 


8 


2006 Sep 17, 18 


21.00 


(10 X 1800 + 1 X 1200) + (6 X 1200 + 1 x 800)s 


191 


157 


1".0 


9 


2006 Sep 18 


21.00 


(12 X 1800)s 


149 


113 


1".0 



Table 4 

Comparison of AAOmega rcdsliift measurements with literature values. The errors presented for Acz and (t(Acz) 
are determined from the 68% confidence limits based on 10000 bootstrap resamplings of the data. 



Literature source 



N, 



match 



Acz 
kms~^ 



(t{Acz) Literature redshift 
uncertainty kms~^ 



Quadrature difference 
kms~^ 



ratio 



Braglia ct al. (2009) 60 

Boschin et al. (2006) 56 

Couch et al. (1998) 29 

Couch & Sharpies (1987) 28 

2df 2003 24 



88 
-18 



+39 
-38 
+23 
24 



-89+41 

,+27 
-26 



-53-^ 



289 
172 
125 
151 
136 



-46 
-40 
19 
18 
+23 
18 
+24 
26 
+20 
23 



276 
90 
90 
100 
50 



86 ± 156 
147 ± 22 

87 ±32 
113 ±32 
127 ±21 



1.4 ±2.6 
2.4 ±0.7 

1.4 ±0.6 

1.5 ±0.7 
2.1 ±0.6 



redshifts where an object has multiple redshift measure- 
ments. Where an object has multiple AAOmega mea- 
surements, we choose the measurement with the highest 
Q value, and where the Q values are equal, we choose the 
measurement with the lower redshift uncertainty. Our fi- 
nal combined catalog contains 1237 extragalactic objects 
with robust redshift measurements within a 15' radius 
of the cluster. The position, rp magnitude, redshift, red- 
shift uncertainty and redshift source for all galactic and 
extragalactic objects in our catalog are tabulated in Ta- 
ble 5. 

3.1.5. Spectroscopic Completeness 

In the forthcoming sections, we will use the redshift 
catalog to search for substructure in Abell 2744. There- 
fore, in order to be confident in the results, it is impor- 
tant that the spectroscopic completeness is well under- 
stood, particularly for the interpretation of maps of the 
projected density of the member galaxies. We use the 
combined redshift catalog to determine the spectroscopic 
completeness, which is defined as the ratio of the num- 
ber of objects with reliable redshift determinations to the 
number of objects in the parent photometric catalog. In 
Figure 8 we present the spectroscopic completeness as a 
function of radius, rp magnitude and a combination of 
radius and rp magnitude. The plots show that we have 
obtained a high level of completeness at all cluster-centric 
radii and for all magnitude ranges. However, the spec- 
troscopic completeness is systematically lower within the 
central ~ 1 Mpc (although we note it is still excellent and 
remains above ^ 70% for all magnitude ranges). This is 
due to a combination of fiber separation limits 30") 
and the high galaxy surface density in the central re- 
gions, which means a small number of objects could not 
be observed. 



Figure 9 shows two images of the spatial distribution 
of the spectroscopic completeness for the magnitude lim- 
its rF< 20.5 (left panel) and rp< 21 (right panel). The 
images were generated by first producing an adaptively 
binned image of the spatial distribution of the galaxies 
in the photometric catalog where each spatial bin con- 
tains ^ 10 galaxies. We made use of the WVT binning 
algorithm by Diehl & Statler (2006), which is a general- 
ization of Cappellari & Copin (2003) 's Voronoi binning 
algorithm. The same spatial binning was then applied 
to the galaxies in the spectroscopic catalog to produce a 
second image. The ratio of the binned spectroscopic to 
photometric images gives the spectroscopic completeness 
map in Figure 9. The maps reveal the high level of spec- 
troscopic completeness we have obtained with our obser- 
vations and, in particular, for rF< 20.5 we obtain > 90% 
completeness over the majority of the 6 Mpc diameter 
field. For rp< 21 we obtain > 80% completeness over 
the majority of the 6 Mpc diameter field. The system- 
atically lower spectral completeness within 1 Mpc noted 
in the radial profile shown in the right panel of Figure 8 
is due to a region of low completeness confined to the 
north west of the center. This patch of lower spectral 
completeness needs to be kept in mind in the forthcom- 
ing sections. 

3.1.6. Cluster Membership 

Robustly distinguishing cluster members from fore- 
ground and background interlopers is a critical first-step 
when conducting a dynamical analysis. Membership al- 
location is a two-step process where the initial step in- 
volves crudely separating cluster members and projected 
interlopers in redshift-space only. This is achieved by 
selecting galaxies within the redshift window of czdus ± 
10000 km s"^ (where Zdus = 0.306; Boschin et al. 2006) 
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Table 5 

Combined redshift catalog. 



R.A. (J2000) 


decl. (J2000) 


rp magnitude 


(cz) 


cz uncertainty 


cz Quality 


C2: source 


deg. 


deg. 




km s ^ 


km 








-00.091094 


20.69 


400971.27 


314.78 


4 


AAi / AAUmega 




-oU.OOo 1 ZZ 


ly.oz 


QQ 1 Q 

yy.io 


1 1 Q QO 


D 


AAI / AAv_/mega 


o.oUoo / 


-oU.4U4oU0 


zU.OO 


c^fim n no 


on nn 





A AT^ / A A r~^m^^.Tr, 

AAi /AAwmcga 


3.304125 


-30.367528 


19.47 


234.52 


65.95 


6 


AAT /AAOmcga 


3.305000 


-30.431861 


20.75 


103.98 


65.95 


6 


AAT/AAOmcga 


3.305750 


-30.351361 


20.48 


28873.99 


20.99 


3 


A AT/ AAOmcga 


3.305833 


-30.339167 


20.27 


92533.37 


23.98 


5 


AAT/AAOmcga 


3.306250 


-30.380945 


22.20 


63794.34 


276.00 


3 


Braglia et al. (2009) 


3.308583 


-30.363083 


19.57 


78171.58 


23.98 


5 


AAT/AAOmcga 


3.310292 


-30.324250 


19.25 


93436.37 


74.95 


4 


AAT/AAOmcga 


3.313542 


-30.398889 


18.91 


91532.63 


23.98 


5 


AAT/AAOmcga 


3.314292 


-30..344611 


18.90 


92477.01 


113.92 


4 


AAT/AAOmcga 


3.313542 


-30.398889 


18.91 


91532.63 


23.98 


5 


AAT/AAOmcga 


3.314292 


-30.344611 


18.90 


92477.01 


113.92 


4 


AAT/AAOmega 


3.314583 


-30.472861 


19.14 


92390.06 


86.94 


4 


AAT/AAOmega 



Note. — This table is available in its entirety in a machine-readable form in the online journal. A portion is shown 
here for guidance regarding its form and content. 



Figure 8. Spectroscopic completeness as a function of magnitude (left) and as a function of radius (right). Also plotted in the right panel 
is the spectroscopic completeness as a function of radius for the magnitude ranges indicated at the lower right of the figure. 



which, as the left panel in Figure 10 shows, defines a 
peak in redshift space associated with Abell 2744. The 
second part of the process refines the membership allo- 
cation by using a slightly modified version of the "shift- 
ing gapper" first employed by Fadda et al. (1996). This 
method utilizes both radial and peculiar velocity infor- 
mation to separate interlopers from members as a func- 
tion of cluster-centric radius (defined with respect to the 
bright cluster galaxy [BCG] closest to the X-ray peak at 
R.A.=00'^14'°20^738, decl.=-30°23'59".90) and is de- 
scribed in detail in Owcrs et al. (2009c). Bricfiy, the 
data arc binned radially such that each bin contains at 
least 40 objects. Within each bin galaxies are sorted by 
their peculiar velocity and interlopers are rejected based 
on the size of the gap in velocity between adjacent (in pe- 
culiar velocity) objects. The limit of 40 objects per bin is 
somewhat subjective and is chosen to minimize both the 
radial annulus used and the likelihood of the algorithm 
rejecting either high velocity bona fide cluster members 
or coherent substructures associated with merger activ- 
ity. Binning by 40 objects produces the best results, as 



shown in the right panel of Figure 10 where we present 
the results of the shifting gapper membership allocation 
procedure. Figure 10 shows that within a cluster-centric 
radius of 3 Mpc the cluster members are well separated 
from the surrounding interlopers. Beyond 3 Mpc the sep- 
aration is less-clcar, so we exclude the galaxies at these 
larger radii during the analysis. After culling the inter- 
lopers, 343 bona fide cluster members remain from which 
we determine a cluster redshift of Zdus = 0.3064 ±0.0004 
and a velocity dispersion of 

^clus — 1497it57kms us- 
ing biweight estimators of location and scale (Beers et al. 
1990), where the quoted errors arc Icr values determined 
using the jackknife resampling technique. 

Also plotted in the right panel of Figure 10 are the 
means (filled circles) and 3a limits (filled diamonds) for 
each radial bin, determined using biweight estimators, 
showing that all excluded interlopers have peculiar ve- 
locities exceeding the Scr limit at the radial bin of inter- 
est. The inadequacy of a simple 3(t clipping using only 
velocity information is highlighted in the right panel of 
Figure 10, where the dashed-lines show that the "icydus 
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limits fail to reject several interlopers at larger cluster- 
centric radii. 

4. SUBSTRUCTURE DETECTION 

The primary motivation for this paper is to detect 
and characterize substructure in order to understand the 
merger history of Abell 2744. While substructure tests 
that use all available information (position on the sky 
and the line of sight velocity) are more decisive (Dressier 
& Shectman 1988; West & Bothun 1990; Pinkney et al. 
1996; Burgett et al. 2004), the results of simpler analy- 
ses are easier to comprehend. In this section, we begin 
with a one-dimensional analysis of the velocity distribu- 
tion, which is most sensitive to detecting substructure 
caused by line-of-sight mergers where the velocity dif- 
ference is large. We then search for substructure in the 
two-dimensional galaxy distribution, which is most sen- 
sitive to detecting well-separated mergers occurring in 
the plane of the sky. Finally, we combine the spatial and 
velocity information to search for local deviations from 
the global velocity distribution which indicate dynamical 
substructures. 

4.1. Substructure in the Velocity Distribution 

Cluster mergers can induce significant distortions in 
the observed velocity distribution which, when a cluster 
is dynamically relaxed, is well approximated by a Gaus- 
sian. Therefore, searching for evidence of substructure 
using velocity information is generally accomplished by 
testing for departures from a Gaussian shape. Here, we 
use the method first outlined by Zabludoff et al. (1993) 
which approximates the velocity distribution by the sum- 
mation of three Gauss-Hermite functions. The method 
is described in detail in Owers et al. (2009c). The co- 
efficient of the first term, /iq « 1, multiplies the best- 
fitting Gaussian contribution to the distribution while 
the second and third terms, with coefficients /13 and /14, 
approximate the third- and fourth-order asymmetric and 



symmetric deviations from a Gaussian shape, similar to 
skewness and kurtosis measurements but less sensitive 
to outliers in the distribution. For Abell 2744, we find 
the coefficients have values of /13 = 0.07 and /14 = 0.06. 
The significance of these values is determined by gener- 
ating 10,000 Monte Carlo realizations of Gaussians with 
N — 343 data points having mean and standard devia- 
tion equal to the best-fitting value derived from the data. 
We find that values of |/i3| > 0.07 occur in only 6% of 
the realizations, while values of |/i4| > 0.06 occur 8% of 
the time. Thus, there appear to be mildly significant 
symmetric and asymmetric deviations from Gaussianity 
present in Abell 2744's velocity distribution. 

The positive value for the term is likely to be due 
to the contribution of galaxies in the cluster outskirts, 
which follow more radial orbits and hence produce a 
more peaked velocity distribution (Merritt 1987), rather 
than the presence of substructure. Because we are inter- 
ested in detecting structure related to a cluster merger 
and not due to orbital anisotropy, we re-determine the 
Gauss-Hermite coefficients using only data from the cen- 
tral 1.5 Mpc region. There are 238 galaxies within this 
region in our sample and for these we measure = 0.12 
and = —0.03. The observed /13 coefficient is highly 
significant with \hs\ > 0.12 occurring in less than 1% of 
the realizations, while the /14 term is not significant and 
I/14I > 0.03 occurs in 47% of the realizations, verifying 
that the mildly significant positive /i4 value measured 
above is indeed due to the contribution of the galaxies at 
larger cluster-centric radii. 

As noted previously, the velocity distribution in 
Abell 2744 is bi-modal (Girardi & Mezzetti 2001; Boschin 
et al. 2006; Braglia et al. 2007). This bi-modality 
causes the significantly positive ft.3 term and Figure 11 
clearly shows a secondary substructure with a peak at 
Upec ~ 2000 km s~^. To further quantify the significance 
of this bi-modality, we use the KMM algorithm (Kaye's 
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Figure 10. Left panel: A histogram of all the reliable redshifts measured within the Abell 2744 field. Arrows mark the initial cut in 
velocity used to crudely isolate the cluster in redshift-space. Right panel: Shift gapper plot for refined member allocation. The black crosses 
represent galaxies allocated as cluster members while the orange open squares are rejected foreground and background galaxies lying close 
to the cluster in redshift space. 
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Figure 11. Left panel: Peculiar velocity histogram for all 343 members. The solid line shows a Gaussian with Vpec = Okm/s and 
(T = 1497 km s~^ and the dashed-line shows the Gauss-Hermite reconstruction of the velocity distribution. Middle panel: Same as the left 
panelfov the 238 members within a l.SMpc radius of the cluster center. Here, the plotted Gaussian (solid line) has a = 1701 hms~^ .Right 
panel: The results of the KMM partition for the members within l.SMpc. The blue Gaussian has Vpec = — 596km/s and a = 1261 kms~^ 
and the red Gaussian has Vpec = 2316km/s and a = 581kms~^. The dashed line shows the combined red and blue Gaussians. 



Mixture Model; Ashman et al. 1994) to fit two Gaus- 
sians to the velocity distribution of the members resid- 
ing within a cluster-centric radius of 1.5 Mpc. As in- 
put for the main component, we set a mean fti^mit — 0, 
dispersion ai^init = 1300 km and estimate the frac- 
tional contribution of the number of galaxies belonging 
to the main component with respect to the entire distri- 
bution is fi^init = 0.7. Likewise, for the secondary com- 
ponent we set fi2,init — 2000 kms~^, <72,init — 500kms~^ 
and f2,init = 0.3. As output, the algorithm returns 
Ml. out = -596kms"\ ai^out = 1261 km s"^ and fi^out = 
0.79 for the main component and ^2,out = 2316 km s~^, 
o'2,out = 581kms~^ and f2,out = 0.21 for the secondary 
component. It is noted that the outputs generated by 
the KMM algorithm are robust to changes in the initial 
input parameters. The resulting Gaussians are overplot- 
ted on the observed velocity distribution in right panel 
of Figure 11. Since we are using the algorithm for a het- 
eroscedastic (i.e., where the variances differ) case, the 
estimate of the significance of the likelihood ratio test 
statistic (LRTS) comparing a bi-modal fit to the null- 
hypothesis uni-modal fit given by the algorithm is unre- 



liable (Ashman et al. 1994). To reliably determine the 
significance of the LRTS, we produce 10,000 randomly 
sampled Gaussian distributions with 238 data points, 
mean -38kms~^ and dispersion 1632 km s~^ (i.e., mean 
and standard deviation which best describes the ob- 
served velocity distribution). For each random dataset, 
two Gaussians are fitted using the KMM algorithm, as 
was done for the observations, where the initial inputs 
are set to the best-fitting outputs for the observed ve- 
locity distribution listed above and the LRTS is com- 
puted. Comparison of the observed LRTS to the dis- 
tribution of 10,000 LRTS produced by the random uni- 
modal Gaussian datasets allows the significance of the 
observed LRTS to be reliably determined. The observed 
LRTS is larger than 99.99% of the LRTSs produced in 
the 10,000 realizations and we conclude that a bi-modal 
fit is significantly better than a uni-modal one. 

4.2. Substructure in the Spatial Distribution of Member 

Galaxies 

Having shown that the velocity distribution in 
Abell 2744 is bi-modal, we now search for further ev- 
idence of merger-related substructure using 2D spatial 
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information. To that end, wc have used the positions of 
the 343 spectroscopicaUy confirmed members to produce 
maps of the smoothed galaxy surface density (Figure 12). 
We used a variable width Gaussian smoothing, where the 
width, a, at each pixel is the projected distance to the 
tenth nearest galaxy, limited to a minimum of 90 kpc and 
a maximum of 500 kpc in high and low density regions, 
respectively. As detected previously (Boschin et al. 2006; 
Braglia et al. 2009), two distinct substructures inhabit 
the central regions of the cluster — the most significant 
one is centered close the cluster center (which was arbi- 
trarily defined in Section 3.1.6 to be the position of the 
BCG nearest the peak in the X-ray surface brightness) 
and a second substructure; ~ 600 kpc to the north. At 
larger scales the distribution of galaxies appears smooth, 
apart from a structure ^ 2.5 Mpc to the west and a mild 
asymmetry to the north-west. 

Given the bimodality detected in the velocity distri- 
bution, we search for any correlation between the two 
velocity peaks and the two substructures seen in the 2D 
galaxy density map for the spectroscopicaUy confirmed 
members (top panels of Figure 12). To achieve this, we 
plot separately the spatial distributions of the galaxies 
allocated to the two KMM velocity partitions in Sec- 
tion 4.1 and smooth the resulting distributions with an 
adaptive circular top hat filter with radius set such that 
the smoothing region contains 10 galaxies. These plots 
are shown in Figure 13 where the left and right panels 
show the distributions of galaxies assigned to the lower 
and higher velocity partition, respectively. The spatial 
distribution of the galaxies assigned to the dominant 
(both in terms of number of galaxies and velocity dis- 
persion) lower velocity partition has two local peaks in 
galaxy surface density. The most significant of these two 
peaks is located coincident with the northern substruc- 
ture seen in Figure 12, while the less significant one is 
approximately coincident with the more central density 
peak seen in Figure 12. The spatial distribution of the 
high velocity partition also has its peak density approx- 
imately coincident with the central density peak seen in 
Figure 12. This indicates that the density peak seen close 
to the center in Figure 12 is not the center of the cluster, 
but is caused by the alignment of two substructures sep- 
arated by a large velocity and projected along our line of 
sight. 

Inspecting the completeness map Figure 9 reveals that 
galaxies in the region of the northern peak are under- 
represented in our spectroscopic sample. Presumably, 
improving the sample completeness in this region would 
enhance the local density peak. Thus the northern peak 
is likely to be even more significant than it appears to be 
in Figure 13. 

The top right panel in Figure 12 shows a close up 
of the smoothed nioinbcr galaxy surface density in the 
central 1x1 Mpc region with Chandra X-ray contours 
overplotted. Neither of the local peaks in galaxy surface 
density coincides with a peak in X-ray surface bright- 
ness. Of particular concern is the lack of any signif- 
icant overdensity coincident with the northwestern in- 
terloper where Boschin et al. (2006) and Braglia et al. 
(2009) both detected substructure using photometrically 
selected cluster members. As noted in Section 3.1.5, 
the region just northwest of the center has the lowest 
spectroscopic completeness and this may be a contribut- 



ing factor in the non-detection of substructure in this 
region. To test this, we obtain the photometric cata- 
log * of Busarello et al. (2002), which covers the central 
2.3 X 2.3 Mpc region, and utilize their publicly available 
photometric redshifts to define cluster membership as 
per the definition in Busarello et al. (2002). The lower 
left panel of Figure 12 shows the adaptivcly smoothed 
galaxy surface density for the photometric members with 
_R < 21 (where 70 < cr < 400 kpc). At this magnitude 
limit, which approximately matches our spectroscopic 
limit of rF=21, there is qualitative agreement with the 
structures detected in the galaxy surface density map 
for the spectroscopicaUy confirmed members. Therefore, 
the non-detection of a significant local overdensity co- 
incident with the northwestern interloper, as detected 
by Boschin et al. (2006) and BragUa et al. (2009), is 
not due to the lower spectroscopic completeness, but 
rather is due to the fact that the majority of the galax- 
ies contributing to the overdensity there are fainter than 
our limiting magnitude. This can be seen in the lower 
right panel of Figure 12 where we have plotted the adap- 
tively smoothed galaxy surface density for the photomet- 
ric members down to the completeness limit of i? = 22.3 
(where 50 < cr < 400 kpc). This map reveals the pres- 
ence of an overdensity approximately coincident with 
the northwestern interloper, consistent with the maps 
of Boschin et al. (2006) and Braglia et al. (2009). This 
result is consistent with Andreon (2001) who, using deep 
if— band imagery, find significantly more dwarf galaxies 
and significantly fewer bright galaxies in this northwest- 
ern region when compared to the central cluster region 
(which is roughly coincident with our defined cluster cen- 
ter). Confirmation that this overdensity is part of the 
Abell 2744 system requires deeper spectroscopic obser- 
vations. 

4.3. Substructure from 3D spatial plus velocity 
information 

We now combine the spatial and velocity information 
to determine if there is any correlation between the sub- 
structures detected in the 2D spatial and ID velocity 
distributions above. This is achieved by making use of 
the K-test which was successfully employed by CoUess & 
Dunn (1996) to detect substructure in the Coma cluster. 
The K-test searches for local departures from the global 
velocity distribution around each cluster member galaxy. 
This is accomplished by using the Kolmogorov-Smirnov 
(KS) test which determines the likelihood that the veloc- 
ity distribution of the n ~ \/N (where N is the number 
of cluster members) nearest neighbors around the galaxy 
of interest and the global cluster velocity distribution, 
which includes the remaining N — n galaxies, are drawn 
from the same parent distribution. The likelihood is de- 
termined by measuring the maximum separation of the 
two cumulative distribution functions. Dobs and deter- 
mining the probability that the Z?-statistic is larger than 
Dobs for the given sample sizes, Pks{D > Dobs)- The 
global measure of the dynamical substructure present 
within the cluster, k, is determined by summing the neg- 
ative log-likelihood for each galaxy over the entire sam- 

* http://vizier.cfa.harvard.edu/viz-bin/VizieR?- 
source= J / A-f A /389/787/ 
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Figure 12. Top left panel: This panel shows the adaptively smoothed galaxy surface density distribution for all spectroscopically confirmed 
cluster members. Top right panel: Close-up of the image in the Top left panel with Chandra isointensity contours overlaid. Crosses mark 
the positions of the 5 bright cluster members. Bottom left panel: Adaptively smoothed galaxy surface density for photometrically defined 
cluster members with magnitudes brighter than R=21 (see text). Bottom right panel: Same as Bottom left panel but for galaxies brighter 
than R=22.3. Note the emergence of a substructure coincident with the X-ray substructure revealed by the Chandra contours. In all 
panels, the color bars show the projected galaxy surface density scale in units of galaxies/Mpc'^. 



pie, i.e. 



(2) 



The significance of k is determined by performing 10,000 
Monte Carlo realizations where, for each realization, the 
observed positions of the galaxies are retained, while the 
galaxy velocities are shuffled randomly, removing any 
correlation between positions and velocities. The ob- 
served K value is then compared to the distribution of 
the 10,000 Monte Carlo realizations of k. A value as 
high as the observed value of k — 481 does not oc- 
cur in the 10,000 realizations, which have a distribution 
that is well modeled by a log-normal distribution with 
/i(lnK) — 4.95 and a{\nK) — 0.19. Thus, the observed k 
value lies 6.6a from the mean of the realizations and 
the upper limit on the probability of observing a value 
this high by chance is lO""*. 
Figure 14 shows the results of the K-test in the form of 



a "bubble plot." At the position of each member galaxy, 
a circle is plotted with radius cx —logPxsiD > Dobs)- 
Clusters of large bubbles reveal regions where the local 
dynamics are different from the global cluster dynamics; 
we define significantly large bubbles as those with radii 
that occur in less than 1% of the 10,000 realizations and 
these are emboldened in Figure 14. The bubbles are also 
color-coded based on the sign of the galaxy's peculiar 
velocity, where red and blue mark positive and negative 
Vpec, respectively. Overplotted are the galaxy surface 
density contours taken from the adaptively smoothed im- 
age shown in the top left panel of Figure 12. There are 
two clear dynamical substructures detected in the central 
1 Mpc; one coincident with the cluster center and an- 
other to the north, coincident with the second projected 
galaxy overdensity. At larger radii, we also detect the 
two substructures to the northwest and south first noted 
by Braglia et al. (2007). 
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Figure 13. Spatial distribution of the two KMM velocity partitions. The positions of the galaxies are marked by open diamonds which 
are color coded according to the KMM partition they were allocated to (Figure 11). The contours show the adaptively smoothed galaxy 
surface density levels, beginning at 5 galaxies/Mpc and incrementing by 10 gals/Mpc. 



Radial gradients in the cluster velocity dispersion are 
expected and, in particular, the velocity dispersion is ex- 
pected to decrease at larger radii due to a combination of 
the radial orbits of infalling galaxies and the decreasing 
escape velocity of the cluster. This can cause spurious re- 
sults in 3D tests, such as the K-test used here, which are 
not related to dynamical substructures (Pinkney et al. 
1996). For this reason, we re-ran the K-test consider- 
ing only those galaxies residing within a central 1.5 Mpc 
radius region. Again, we find the observed value of 
K = 200 does not occur in any of the 10,000 Monte 
Carlo realizations and lies ^ 3.6cr from the mean of the 
simulated k distribution, which has //(Iuk) = 4.58 and 
cr(lnK) = 0.20. 

The K-test is an excellent tool for locating dynami- 
cally distinct substructure, however it provides no insight 
into the origin of the differences in the kinematics of the 
galaxies within the substructures. As a first step towards 
visualizing and characterizing the substructure revealed 
by the K-test in Figure 14, we present Figure 15 which 
shows the line-of-sight velocity field. The value of each 
pixel in the line-of-sight velocity field map is the trimmed 
mean of the velocity distribution for the 15 nearest neigh- 
bors in projection. Only pixels where the radius to the 
15th nearest neighbor is < 1000 kpc are presented in Fig- 
ure 15. The trimmed mean is measured by first deter- 
mining the median and median absolute deviation (which 
are insensitive to outliers) as proxies for the mean and 
standard deviation for the distribution, trimming those 
points which are further than Bcr from the median and 
calculating the mean of the remaining points. The ve- 
locity field is complex and, most significantly, there is 
clear evidence for velocity structure within the central 
IMpc region where a low velocity (~ — ISOOkms"^) 
structure is seen to the north, while a high velocity struc- 
ture (~ 2400 km s~^) is seen to the south. These two 
structures coincide with the two substructures revealed 
by the K-test in Figure 14. 



We now attempt to disentangle these substructures 
and ascertain their dynamics. To do this, we again uti- 
lize the KMM algorithm of Ashman et al. (1994) as in 
Section 4.1, however, here we use both spatial and veloc- 
ity information to partition the cluster into its dynami- 
cal substructures. Since we are mainly interested in the 
dynamics of the central region, where we have Chandra 
data, and the KMM algorithm can be biased by outliers, 
we consider only members within a cluster-centric radius 
of 1.5 Mpc for the KMM analysis. As input, the KMM 
algorithm requires the number of partitions needed to 
describe the data, as well as the approximate mixing 
proportions and estimates of the means and covariance 
matrices for each partition. One of the major uncer- 
tainties when using this method is in determining the 
number of partitions required to adequately describe the 
data. We implement two strategies to overcome this un- 
certainty. First, we use the results from the re-test, in 
combination with the galaxy surface density maps, to 
guide us in selecting the number of likely substructures. 
Second, we utilize the hypothesis testing capabilities of 
the KMM algorithm by measuring the LRTS, which gives 
a measure of the improvement in the fit for a complex, 
g-partition, model over a simpler, go-partition model. In 
Ashman et al. (1994), the simpler model was a unimodal 
one (i.e., go = 1) where the significance of the LRTS, ex- 
pressed as a P- value, was estimated by comparison with 
a distribution, which is only suitable when compar- 
ing two homoscedastic (i.e., the same variance), univari- 
ate models — clearly not the case here. Furthermore, we 
would like to assess how significantly the fit is improved 
by the addition of one partition, i.e., the f;-partition case 
over the go — g — 1-partition case, in order to decide how 
many partitions to use. Therefore, we employ a para- 
metric bootstrap resampling technique to determine the 
P-values as follows. We use KMM to fit the go- and g- 
partition 3D Gaussian models to the data. We then pro- 
duce 5000 bootstrapped samples by randomly sampling 




Figure 14. Bubble plot outputs from the k test. The bold bubbles are those deemed to be significant insomuch as they only occur in 1% 
of 10,000 realizations. The total k test is significant at the 6.6a- level and the observed value does not occur in all 10,000 realizations. Blue 
and red bubbles have negative and positive Vpec, respectively. The contours are galaxy density contours generated from the adaptively 
smoothed image in the top left panel of Figure 12 and are linearly spaced by 10 in the interval 5-245 gals Mpc~^ . Note the two clusterings of 
significant bubbles are coincident with overdensities in the projected galaxy density. The inset panels show the peculiar velocity histograms 
for the regions indicated by the green circles, which have radii of 350 kpc. The arrows mark the velocities of the bright ('rp< 17.8) member 
galaxies within each green region. 



the best fitting go-partition 3D Gaussian mixture model. 
To each of these bootstrap samples, we re-fit the go- and 
g-partition models, using the means and covariances of 
the models fitted to the data as initial input estimates 
for the KMM algorithm and remeasure the LRTS. The 
observed value for the LRTS is then compared to the null 
distribution of bootstrapped LRTS values to determine 
the P-value. 

Figure 14 reveals two clear substructures: one with 
Vpec ^ 2300 kms~^ close to the cluster center, which we 
label Clump A and one with Vpec ^ —1600 kms~^ and 



^ 700 kpc north of the cluster center, which we label 
Clump B. There is also a possible substructure approx- 
imately spatially coincident with Clump A, but with a 
Vppc similar to that of Clump B, which we call Clump C. 
These partitions were used as input for the KMM al- 
gorithm. Thus, we fit three models: a = 2-, 5 = 3- 
and g — 4-partition model and at each stage we deter- 
mine the P-value for the null hypothesis go =5 — 1 
fit, as outlined above. To determine the initial esti- 
mates for input into the KMM algorithm, we take ad- 
vantage of the fact that each of these clumps houses a 
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Figure 15. The line-of-sight velocity field determined at each 
pixel using the trimmed mean of the velocity distribution of the 
15 nearest galaxies. The black contours show the galaxy density 
isopleths and the pixel size is 50kpc. 

bright {rp< 17.8) galaxy which has Vpec at the peak of 
the velocity distribution (Figure 14). We assume that 
these bright galaxies lie at the spatial and dynamical 
centers of each partition and crudely define partition 
membership as any galaxy within Vpec = ±1500 kms~^ 
and radius < 350 kpc of the central galaxy. The means 
and covariance arrays for these initial partitions serve 
as the initial estimates for the KMM algorithm. The 
results of the g = 2,3 and 4-partition KMM fits are pre- 
sented in Figure 16. The P — values for these fits are 
Pg=2 < 0.0002, Pg^3 = 0.02, and Pg=i = 0.93. Thus, 
the 5 = 3 partition fit, with Clumps B and C combined, 
is the statistically preferred one, despite a spatial sepa- 
ration of ^ 700 kpc between the two bright galaxies in 
these clumps. 

5. THE NATURE OF THE DETECTED SUBSTRUCTURES. 

In this section, we discuss the relationship between the 
structures detected in the X-ray and optical analyses and 
present interpretations of the nature of these structures. 
As a visual aid to the discussion below, we present Fig- 
ure 17 which is an i?— band image taken with the AAT. 
Overlaid are contours from the X-ray residual map (Fig- 
ure 2) showing the relevant X-ray substructures. The 
cluster member positions are overplotted and color coded 
to match the allocations from the KMM g = i partition 
(Section 4.3). While this partition was not found to be 
favored over a g = 3 partition fit, it serves the purpose 
of a visual aid by highlighting the dynamical structures 
discussed below in a more coherent manner. 

5.1. The Northern core and northern major remnant 

core 

In their analysis of the dynamics of Abell 2744, Boschin 
et al. (2006) detect the northern major remnant core 
in both the 2D and 3D analyses but discard the hy- 
pothesis that it is the undisturbed remnant core of the 



main cluster. The combination of the new Chandra X- 
ray data and comprehensive AAOmega spectroscopy re- 
veal that the substructure north of the central regions of 
Abell 2744 is more significant than previously thought. 
The strongest evidence for this assertion comes from Fig- 
ure 13 which shows that when the contribution to the 
galaxy surface density from the high velocity southern 
component is removed, the most significant peak in the 
galaxy surface density is coincident with the northern 
substructure. In Section 4.3, we show that this galaxy 
surface density enhancement is detected as a dynami- 
cally distinct substructure with Vpec ^ —1600 km s~^ and 
ay ~ 800kms~^. Furthermore, the deeper Chandra ob- 
servations resolve a gas substructure and reveal a tail of 
gas trailing to the south which curves towards the west 
(Figure 1). The abundance map presented in Section 2.3 
(Figure 4) shows that this gas substructure has signif- 
icantly higher metallicity than the global cluster value, 
while the analysis of detailed spectra (NCI in Table 2) 
show that this gas structure has a significantly lower tem- 
perature {kT ^ 7keV) than the global cluster tempera- 
ture. 

Observations of relaxed clusters show that the central 
core regions are characterized by a steep decline in gas 
temperature and a corresponding increase in metallicity 
(Vikhlinin et al. 2005). Given these observations and 
the thermodynamic properties outlined above, we inter- 
pret the northern gas substructure as the remnant of a 
cool-core which has survived the merging activity. In 
addition to this, the large velocity dispersion, the signif- 
icant number of galaxies in this region, along with the 
fact that this structure harbors one of the bright clus- 
ter members (with Vpec ^ — 1200 km s""*^), we postulate 
that this northern component has a significant mass and 
is the core of the more massive system involved in the 
central merger in Abell 2744. 

5.2. The central tidal debris 

This structure, located coincident in projection with 
the high velocity southern substructure (Section 5.3), but 
separated in peculiar velocity by ~ 4000 km s""'^ as dis- 
cussed in Section 4.3 (Figure 14 inset), has previously 
been interpreted as being the core of the main negative 
Vpec component in the Abell 2744 system (Kempner & 
David 2004; Boschin et al. 2006). This interpretation 
is supported by the local peak in the surface density 
of the galaxies and the presence of a luminous galaxy 
with Vpec ~ — 1600kms~^. However, as discussed in Sec- 
tion 4.2 (see also Figure 13), the peak in the projected 
galaxy density here is artificially enhanced by the high 
velocity southern substructure. When this is accounted 
for, it can be seen that the more significant peak in the 
galaxy surface density of the negative Vpec component 
is in fact the northern substructure, which we interpret 
as the remnant core of the main cluster, as discussed in 
Section 5.1. 

Despite not being the core of the negative Vpec compo- 
nent, this structure is a dynamically distinct component 
with a Vpec ~ — 1600kms~^. In Section 4.3, we found 
that adding a component representing this structure to 
the KMM analysis failed to improve the fit significantly 
over the simpler 3-partition fit. Therefore, we tentatively 
identify this structure as tidal debris stripped from the 




Figure 16. These plots show the results of the KMM partitioning 



of the data using the full spatial plus velocity information. 
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Figure 17. This image shows an AAT R— band image of Abell 2744. Contours from the X-ray residual map (Figure 2) are shown in 
black. The regions show cluster members and are color/shape coded with red diamonds, green diamonds, blue squares and cyan circles 
corresponding to members allocated to KMM partitions labeled in the lower left panel of Figure 16 as Clump B (Northern major remnant 
core), Clump C (Central tidal debris), Clump A (Southern minor remnant core) and the Main Cluster. 



23 



main cluster (the remnant core of which now has mainly 
to the north) during the core passage phase of the merger 
with the high velocity southern substructure. 

5.3. The Southern compact core and southern minor 
remnant core: A Bullet-like remnant viewed from 
nearer to the merger axis 

The southern compact core stands out as a prominent 
surface brightness enhancement in the Chandra images 

(Figures 1 and 2). It is conspicuous in the thermody- 
namic maps presented in Section 2.3 as a region with 
significantly cooler, higher metallicity and lower entropy 
gas (Figures 3 and 4). This provides strong evidence for 
the southern compact core being the remnant cool core 
of a merging substructure. In Section 4.1 we found that 
the cluster peculiar velocity distribution was bimodal, 
containing a high velocity component. In Section 4.2 
we explored the spatial distribution of this high velocity 
component, finding that it is spatially compact, while 
the combination of spatial and velocity information in 
Section 4.3 confirmed the existence of this structure as a 
dynamically distinct, spatially compact entity which we 
label as the southern minor remnant core in Figure 17. 
The proximity of the southern minor remnant core to the 
southern compact core indicates that they are associated; 
the bright galaxy accompanying the southern minor rem- 
nant core is located only ^ 140 kpc to the southeast of 
the X-ray peak due to the southern compact core. The 
peculiar velocity of the southern minor remnant core, ob- 
tained with the KMM analyses in Sections 4.1 (for the 
ID velocity analysis) and 4.3 (for a full 3D velocity plus 
spatial information analysis), is Vpec = 2300 — 2500km/s. 
We find evidence that the southern compact core also 
has a high peculiar velocity by considering the shock 
front to the southeast (Section 2.4.2). The Rankine- 
Hugoniot shock jump conditions can be used to derive 
a Mach number for the shock from the density jump and 
thus a shock velocity (e.g., Eqn 14 from Markevitch & 
Vikhlinin 2007; Landau & Lifshitz 1959). For the mea- 
sured density jump of l.SOlo;^!, we determine a Mach 
number of M = 1.41+Q jg^. We can check this result 
by using Eqns 4 and 14 from Markevitch & Vikhlinin 
(2007) and the temperature measurements taken on ei- 
ther side of the front to obtain M = I.8IIQ37, consistent 
within the errors. The sound speed for a ~ 8.6 keV gas is 
^ 1523 km/s, therefore the shock velocity is ~ 2150km/s, 
using the Mach number derived from the density jump. 

We note that this method does not give the total ve- 
locity of the shock, rather it gives the shock velocity in 
the plane of the sky. This is due to the following. In 
order that the shock front be observable as an edge in 
the surface brightness, our line of sight miist be tan- 
gent to some point on the shock front surface. Assuming 
we know the radius of curvature, we measure the den- 
sity jump at the point which is tangent to our line of 
sight. The Rankine-Hugoniot shock jump conditions re- 
late the measured density jump to the Mach number of 
the component of gas velocity travelling perpendicular 
to the shock surface. Since our line of sight is tangent 
to the point at which we measure this Mach number, 
the measured velocity component must be in the plane 
of the sky. We can use this fact, in combination with 
the line of sight velocities of the northern major rem- 



nant core and the southern minor remnant core, to con- 
strain the inclination of the merger axis along the line 
of sight and to approximate the shock velocity. To do 
this, we assume that the southern compact and north- 
ern cores have the same line of sight velocities as their 
respective galaxy components, i.e., the southern minor 
and northern major remnant cores (measured from the 
3D KMM partition in Section 4.3). The shock leading 
the southern compact core is meeting the unshocked gas 
associated with the northern core with a line of sight ve- 
locity V108 — 2574-1-1658 = 4232 km s~^ which, assuming 
the unshocked gas has kT^ 8.6 keV, gives a line of sight 
Mach number of Mios ^ 4232/1523 = 2.78. The inchna- 
tion angle of the velocity vector to our line of sight is then 
arctan(1.41/2.78) ~ 27° and the total Mach number is 
Mtot = 3.1 leading to a shock speed of ~ 4750 kms~^, 
similar to that observed in the Bullet cluster (Markevitch 
et al. 2002; Markevitch 2006). Using the Mach number 
derived from the temperature jump gives a similar value 
for the inclination angle of ^ 33°, the total Mach num- 
ber, Mtot — 3.31 and a shock speed of 5041 kms^^. We 
note that, as in the Bullet cluster (Springel & Farrar 
2007; Mastropietro & Burkert 2008), the shock velocity 
is measured relative to the surrounding gas which will 
have high bulk motions due to the merger. These mea- 
surements also rely critically on the assumption that the 
velocities of the gaseous and galactic components is the 
same. Future X-ray telescopes, such as the International 
X-ray Observatory, will have high enough spectral reso- 
lution to directly probe the line of sight gas kinematics 
and will enable us to test this assumption. 

The measured velocity dispersion of the southern mi- 
nor remnant core ranges from 581 km/s, measured from 
the ID KMM partition, to 441 km/s, measured from 
the 3D KMM partition. This dispersion is likely to be 
an underestimate of the real dispersion due the KMM 
partitioning which allocates objects which are probably 
physically associated with the southern minor remnant 
core, but in the tails of the velocity distribution, to the 
main cluster. This eff'ect is seen in the right-most panel 
of Figure 11, which shows that objects with peculiar ve- 
locities overlapping with the main structure have been 
allocated to the main structure, thereby truncating the 
southern minor remnant core's velocity distribution. The 
act of merging will also aff'ect the measured velocity dis- 
persion, since the merger acts to disperse the southern 
minor remnant core's galaxies (Pinkney et al. 1996). The 
galaxies most likely to be stripped are those that are least 
tightly bound to the southern minor remnant core, e.g., 
those in the outskirts, or those in the tails of the velocity 
distribution. This stripping, compounded by the lack of 
certainty in the allocation of galaxies as members of the 
southern minor remnant core, mean it is extremely dif- 
ficult to gauge the southern minor remnant core's initial 
velocity dispersion. In any case, the measured veloc- 
ity dispersion indicates the southern minor remnant core 
contains significant mass, as does the measured temper- 
ature of ~ 7keV for the southern compact core. 

Given the above evidence, we interpret the southern 
compact core as a Bullet-like system (Markevitch et al. 
2002) merging along an axis that lies well out of the plane 
of the sky. To support this interpretation, we present 
Figure 18, which shows a model for the X-ray appear- 
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ance of the Bullet cluster (as used in Bradac et al. 2006) 
viewed from a range of directions. Briefly, the Chandra 
X-ray image of the Bullet cluster is fitted with a density 
model which consists of a shuttlecock with concave sides 
(for the Bullet substructure), power-law profiles across 
the cold and shock fronts, a flattened pancake (for the 
main cluster's emission) and a beta-model at larger radii. 
The XI M software (Heinz & Briiggen 2009) is used in 
combination with this density model to simulate 0.6-5 
keV Chandra images with 100 ks exposures and projected 
so that the merger axis is inclined to our line of sight with 
angles of 0i„c = 90, 15, 35, 48, and 60. To best match the 
orientation of the southern compact core in Abell 2744, 
the images shown in the lower 4 panels of Figure 18 are 
again rotated in the plane of the sky. 

Comparing the simulated images in Figure 18 with the 
X-ray morphology of the southern compact core, we sur- 
mise that an inclination angle in the range 9inc = 35—60° 
provides the best qualitative fit. This is based primarily 
on the observability of the shock front (Section 2.4.2) and 
on the blunted cone appearance of the southern compact 
core in the Chandra images (Figure 1). These inclina- 
tion angles compare well with that derived above from 
the combination of the shock and kinematic parameters, 
particularly when considering uncertainties inherent in 
both techniques. 

5.4. The northwestern interloper 

In their analysis of the northwestern interloper, KD04 
presented tentative evidence for a cold front and a bow 

shock on the eastern side of the northwestern interloper 
(their Figure 3.). In our deeper observations, we find no 
convincing evidence for the existence of a shock front at 
the position and orientation indicated in KD04. How- 
ever, both of the putative fronts appear to be part of a 
larger front on the northeastern side of the northwestern 
interloper, as indicated in Figure 1. The analysis in Sec- 
tion 2.4.3 reveals that this is a cold front. In Figure 2 
we noted that the emission surrounding the northwest- 
ern interloper extends towards the south. The analyses 
conducted in Sections 2.3 and 2.4 indicate that this ex- 
tension is cooler and has lower entropy than the ICM 
immediately to the east and west. We interpret the low 
entropy gas as a trail of gas stripped from the northwest- 
ern interloper by ram pressure. 

While we detect no dynamical system associated with 
the northwestern interloper, the pair of bright galaxies ly- 
ing ^ 140 kpc and ^ 320 kpc southeast of the X-ray peak 
have peculiar velocities of 561 kms^^ and 663 km s^^, re- 
spectively. The velocity field (Figure 15) has values of 
400kms~^ in the region of the northwestern inter- 
loper, consistent with the velocities of the pair of bright 
galaxies. This hints at a dynamically distinct system. 
However, because the velocities are close to those of the 
main cluster and because of the small number of spec- 
troscopically confirmed cluster members in this region, 
dynamical substructure is more difficult to detect at a 
significant level (Pinkney et al. 1996). We note also that 
if the northwestern interloper has passed pericenter, as 
we argue in Section 6.2, it is plausible that a large frac- 
tion of the galaxies have been tidally stripped. This will 
further inhibit the detection of dynamical substructure 
associated with the northwestern interloper. 

Using the photometric redshift catalog of Busarello 



et al. (2002) to define cluster members, probing ^1.3 
mag fainter than our spectroscopy, we do detect a local 
peak in the galaxy surface density coincident with the 
northwestern interloper's X-ray peak. We showed that 
the non-detection of a local peak in the 2D galaxy distri- 
bution of the spectroscopically conflrmed cluster mem- 
bers is not due to spectroscopic incompleteness, but due 
to the majority of galaxies associated with the overden- 
sity detected using the photometrically defined members 
having magnitudes fainter than our limiting magnitude 
for the spectroscopy. Thus, spectroscopic observations 
probing to fainter magnitudes are required to confirm 
the existence of a substructure associated with the north- 
western interloper. 

6. A MERGER SCENARIO 

Here we present a merger scenario based on our new 
observations. Our interpretation of the data reveal a 
merger scenario which agrees with previous authors in 
the broadest sense, i.e., we are observing a major merger 
in the central regions of the cluster and a more minor 
infall of a group to the northwest. However, the details of 
our merger scenario differ from those outlined previously 
in KD04 and Boschin et al. (2006). In the following, 
we sketch an outline for this major-plus-minor-mergcr 
scenario, but stress that confirmation of our hypotheses 
require rigorous testing in the form of planned Abell 2744 
specific merger simulations. 

6.1. The central major merger 

The central region of the cluster has extremely complex 
X-ray morphology, temperature and metallicity struc- 
ture, as well as multiple dynamically distinct compo- 
nents. These features can be reconciled with an off cen- 
ter, two body merger after core passage, with the merger 
axis nearer to our line of sight than the plane of the sky. 
Below, we provide evidence for these assertions. 

During cluster mergers, the strong ram pressure felt 
by the gas means that it can become decoupled from the 
coUisionless galaxies which will lead the gas (e.g., Bar- 
rena et al. 2002; Markevitch et al. 2002; Clowe et al. 
2006; Mahdavi et al. 2007; Bradac et al. 2008). From 
the observed offset between the galaxies and gas, the 
X-ray morphologies, and also the velocities of the sub- 
clusters, we can determine the direction of motion of the 
two subclusters. The northern major remnant core is 
blueshifted with respect to the systemic velocity of the 
cluster complex (fpec — — 1600kms~^), the galaxy sub- 
structure is offset from the northern core to the north- 
west, with the bright cluster galaxy associated with the 
remnant core leading by ~ 100 kpc. This evidence indi- 
cates that some component of the northern subcluster's 
motion is directed towards the north- northwest, while a 
large portion is also directed towards us along the line 
of sight. Further evidence for the north- northwest mo- 
tion is revealed by the enhanced X-ray emission to the 
southeast of the northern core (Figures 1 and 17), which 
we interpret as a trail of stripped gas. The curvature of 
the gas trail implies significant westward acceleration, so 
that the northward moving core must have passed to the 
east of the core of the other major subcluster, which we 
associate with the southern compact core. 

Using a similar argument, we can ascertain the direc- 
tion of motion of the southern subcluster. The brightest 
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Figure 18. Top left: A2744 Chandra image. Top Right: Simulated lOOks Chandra ACIS-I observation of tiie Bullet cluster gas model (at 
z=0.308) viewed in the same projection as the Bullet cluster, i.e., 9i„c = 90 where 9i„c gives the angle of the merger axis with respect to 
the line of sight. Middle left: As for top right with 9inc = 15. Middle Right: As for top right with 9i„c = 35. Bottom Left: As for top right 
with 9i„c = 48. Bottom Right: As for top right with 9inc = 60. The bottom four panels are rotated in the plane of the sky to best match 
the orientation of Abell 2744 shown in the top left panel. 



galaxy associated with the southern minor remnant core 
is offset some ^ 140 kpc south of due east of the peak 
in the X-ray emission associated with the southern com- 
pact core and assuming that this bright galaxy was once 
at the center of the subcluster and coincident with the 
peak in the X-ray emission, this offset indicates some 
fraction of the southern subcluster's motion is directed 
towards the south-east. Further evidence for motion to- 
wards the southeast comes from the Chandra image in 



the form of the edge to the southeast (Figure 1), which is 
part of the shock front (Section 5.3), and the morphology 
of the southern compact core, which exhibits a blunted 
cone structure where the compact core is the head of 
the blunted cone, and the tail of stripped gas to the 
northwest forms a wedge-like structure. The mean pecu- 
liar velocity of the southern minor remnant core galaxies 
(^ 2300 — 2500 km s~^) indicates that a large fraction of 
the southern substructure's motion is directed away from 
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us along our line of sight. However, the edge and blunted 
cone morphology of the southern compact core show that 
the direction of motion cannot be entirely directed along 
our line-of-sight. 

These arguments require that the northern and south- 
ern subclusters have passed one another and now are 
moving apart. The gaseous envelopes of the two subclus- 
ters have largely been stopped and stripped from their 
former hosts by shocks, separating them from the rem- 
nant cool cores, which can survive the off-center merger 
due to their higher densities and pressures. In our merger 
scenario, the central peak of emission is the gas which 
has been stripped from the outer parts of the two sub- 
clusters (luring tlw; core passage phase of the merger. 
This interpretation is supported by high temperatures in 
the central gas, consistent with shock heating and adi- 
abatic compression during the core passage. Patches of 
low metallicity gas suggest that at least some of this gas 
originates from outside the more highly enriched cores 
(Vikhlinin et al. 2005), as described above. 

According to the KMM partitioning in Section 4.3, a 
majority (188 out of 238) of the galaxies are allocated to 
a partition with jj, = — 105kms~^ and a = 1438 kms~^. 
This is also clearly seen in Figure 10 where within a 
cluster-centric radius of 500 kpc there is a clear seg- 
regation in velocity with essentially no galaxies with 
Vpec = Okms~^, whereas at i? > 500 kpc the veloc- 
ity distribution appears to be unimodal and centered 
around Wpec ~ Okms""'^. Presumably the galaxies at 
larger radii were once part of the less-tightly-bound out- 
skirts of the southern and northern substructures. The 
outskirts of the two merging clusters decoupled from the 
tightly bound cores during pericenter and are in the pro- 
cess of mixing to form the outskirts of the new cluster. 

In our interpretation, the northern subcluster was the 
more massive of the two merging systems and the soiith- 
ern subcluster is the remnant of the less massive one, con- 
sistent with the higher velocity dispersion of the north- 
ern major remnant core. However, the southern compact 
core is brighter now. This suggests that the smaller sub- 
cluster hosted a more robust cool core, i.e. a greater 
mass of low entropy gas, prior to the merger. The higher 
entropy northern core of the larger subcluster would then 
be more prone to disruption during the merger. Indeed, 
the clumpy tail of gas trailing the northern core indicates 
that a significant fraction of its gas has been removed and 
it is in the process of being destroyed by the merger. 

6.2. The infall of the northwestern interloper 

The deeper Chandra observations presented here fail to 
confirm the proposal by KD04 that the northwestern in- 
terloper is traveling eastward on its first passage through 

the cluster. KD04 used the positions of two bright clus- 
ter galaxies ~ 140kpc and ^ 320kpc southeast of the 
northwestern interloper, along with tentative detections 
of cold and shock fronts on the eastern side of the inter- 
loper, to support this scenario. We find little evidence 
for the shock proposed to exist by KD04, while their pro- 
posed cold front is part of a more extended cold front. 
Furthermore, it is difficult to explain the large separation 
of gas and galaxies prior to a core passage in a siibchister 
which appears to have both significant mass (indicated 
by its temperature kT~ 5keV) and ICM density (indi- 
cated by its X-ray brightness). 



Three main features in the deeper X-ray data provide 
support for an alternative interpretation. They are 1) 
the cold front to the north-east, 2) the extended region 
of enhanced X-ray surface brightness to the south and 
3) the swirl of low entropy gas seen in Figure 3, which 
is associated with the southern extension and is curved 
in a southeasterly direction towards the cluster center. 
The cold front to the north-east, along with the exten- 
sion of the emission towards the south, indicate that the 
northwestern interloper is currently traveling towards the 
north-northeast. The curvature of the swirl of low en- 
tropy gas indicates significant angular momentum. Us- 
ing these features to guide our interpretation, we suggest 
that the northwestern interloper fell into the main cluster 
from the south or southeast, initially traveling roughly 
north or northwest and passing the main cluster core off- 
center to the southwest. The gravity of the main clus- 
ter has subsequently deflected it onto its current path 
towards the north- northeast (see, e.g., the evolution of 
the 3:1 and 10:1 mass ratio, offcenter mergers in Poole 
et al. 2006). The fact these features are observed at all, 
along with the non-detection of any significant dynam- 
ical substructure in this region, indicates the merger of 
the northwestern interloper is taking place fairly close to 
the plane of the sky. We note, however, that the non- 
detection of dynamical activity here may be due to the 
limited depth and completeness of our sample. 

One issue with the above interpretation for the north- 
western interloper is the position of the two bright galax- 
ies which are offset to the southeast from the peak in the 
X-ray emission. Assuming these bright galaxies trace 
the dominant mass component (i.e., the coUisionless dark 
matter core) of the northwestern subcluster, we would 
expect to see the classic Bullet-like scenario where the 
coUisionless dark matter and galaxies lead the coUisional 
gas component and thus to see the bright galaxies offset 
to the northeast of the gas clump. However, we note that 
during the later stages of a merger, the simulations of 
Mathis et al. (2005) showed that the cool, dense cores of 
subclusters can lead the dark matter cores. Indeed, this 
phenomenon has been observed in the clusters Abell 168 
(HaUman & Markevitch 2004) and Abell 754 (Marke- 
vitch et al. 2003). The physical mechanism producing 
this offset is the rapid changes in ram pressure felt by 
the subcluster gas as it traverses the main cluster; dur- 
ing closest approach, the ram pressure peaks and the gas 
and dark matter separate, with the dark matter leading. 
After the core passage phase, the ram pressure drops 
significantly due to the decrease in the surrounding ICM 
density and the cool core begins to fall back towards its 
parent dark matter core, eventually overshooting it in 
a "ram pressure slingshot" as explained in Markevitch 
& Vikhlinin (2007). This effect is also simulated in As- 
casibar & Markevitch (2006) and their Figure 12 shows 
qualitative agreement with what is proposed above and 
with the temperature morphology of the northwestern 
substructure. High resolution weak lensing maps are re- 
quired to determine the location of the local peak in the 
mass density corresponding to the northwestern inter- 
loper and to test the scenario outlined above. 

While the northwestern subcluster is due to a more mi- 
nor merger than the central merger, there should still be 
observable effects due to the core passage, particularly in 
the ICM. Although the central major merger complicates 
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this matter and may erase any effects due to the north- 
western sul^chister, we tentatively point out the bridge of 
emission in Figure 1 along with the slight asymmetry of 
the central cluster emission in the direction of the north- 
western interloper, as evidence of the previous interac- 
tion with the central ICM. We also suggest that the hot 
regions lying between the main cluster and the north- 
western interloper (Figure 3) were produced by shocks 
and/or compression of the ICM during the subclusters 
passage. 

6.3. Comparison to previous interpretations 

The interpretations outlined above for the different 
structures revealed by our new observations give cause to 
re-visit the hypotheses outlined in Boschin et al. (2006) 
and KD04 for the merging history of Abell 2744. Con- 
sidering the central regions, KD04 propose a merger sce- 
nario where Abell 2744 is seen just prior to a core passage 
and the merger is occurring along a north-south direction 
with the majority of the motion directed along our line- 
of-sight. Boschin et al. (2006), using a slightly larger 
sample of spectroscopically confirmed cluster members, 
agreed with KD04 that the central merger is occurring 
along a north-south direction with the majority of the 
motion directed along our linc-of-sight, but suggested 
that the merger is in a more advanced post-core-passage 
phase. Both KD04 and Boschin et al. (2006) suggest that 
the northwestern subcluster is infalling onto the main 
cluster from the west and can be treated as separate from 
the major merger occurring in the central regions. 

With regard to the southern subcluster, our interpre- 
tation of the new observations is mostly in accordance 
with the conclusion of Boschin et al. (2006) that it is 
viewed after a core passage, has a large fraction of its 
motion directed along the line-of-sight and some compo- 
nent to the south. We slightly modify this scenario based 
on the detection of the shock front and the blunted cone 
X-ray morphology, which indicate a more southeasterly 
direction of motion. The detection of the shock front also 
indicates that the component of the velocity perpendic- 
ular to our line-of-sight must be larger than previously 
suggested. Indeed, our observations of the shock front, 
combined with the dynamical analyses, allow us to con- 
strain the direction of motion of the southern subcluster 
to be inclined at around 30° to the line of sight. 

Our scenario differs most from that of Boschin et al. 
(2006) in the importance of the northern subcluster. 
Boschin et al. (2006) interpret this structure as an addi- 
tional feature in their scenario due to the lack of a coun- 
terpart in the X-ray images. Our deeper Chandra data 
reveals the northern core as a cool, high metallicity struc- 
ture which is trailing the dynamically distinct northern 
major remnant core galaxies. Furthermore, the larger 
sample of spectroscopically confirmed cluster members 
provide strong evidence that, when separated by veloc- 
ity, the most significant projected galaxy overdensity is 
the northern major remnant core. Taking these new re- 
sults in concert provides strong evidence in support of our 
view that the northern subcluster is the surviving core 
of the main component involved in the central merger. 
Thus, while we agree with Boschin et al. (2006) that 
we are observing a post-core passage merger in the cen- 
tral regions, we suggest a somewhat later phase, given 
the large (^ 800 kpc) projected separation of the main 



northern subcluster from the high velocity southern sub- 
cluster. 

Turning to the interloping northwestern subcluster, the 
deeper Chandra data have muddied the waters somewhat 
in terms of a simple interpretation and indicate that the 
direction of motion of this subcluster is currently to the 
north/northeast perpendicular to the direction of mo- 
tion suggested by both KD04 and Boschin et al. (2006). 
Both Boschin et al. (2006) and Braglia et al. (2009) de- 
tect an overdensity of red galaxies associated with the 
bright cluster galaxy lying ~ 140 kpc to the southeast 
of the X-ray peak. We propose that this offset has oc- 
curred due to a ram pressure slingshot which caused the 
gas to "overshoot" the galaxies after pericentric passage. 
This interpretation provides a better fit to the X-ray data 
and we note that the galaxy overdensities found here and 
in the Boschin et al. (2006) and BragUa et al. (2009) 
studies rely on photometric cluster member allocation, 
rather than the more robust spectroscopic identification 
of cluster members, which is less prone to projection ef- 
fects. A more rigorous investigation involving deeper 
spectroscopy and detailed simulations are required to test 
the different infall scenarios for this substructure. 

7. SUMMARY AND CONCLUSIONS 

We have presented an analysis of the merging chis- 
ter Abell 2744 based on new Chandra X-ray data and 
AAOmega optical spectroscopy. The deeper Chandra 
data reveal a plethora of substructure including: 

1. A significant structure to the north with a trail of 
X-ray emission to the south and curving westward. 
This northern core harbors cooler gas with higher 
metallicity than its surroundings, indicating it is a 
remnant cool core. 

2. The Bullet-like southern compact core which is sur- 
rounded by hot, shocked gas. There is an edge to 
the southeast of the southern compact core, which 
has the characteristics of a shock front with a Mach 
number M = lAlto.H- 

3. A global peak in the X-ray emission with no corre- 
sponding galaxy overdensity, which is interpreted 
as the stripped atmospheres of the merging sub- 
clusters. 

4. The northwestern interloper has an edge to the 
northeast and a tail of emission extending towards 
the south. Analysis of the spectra of these features 
indicate that the edge is a cold front and that the 
tail of gas harbors cooler, lower entropy gas when 
compared to the surrounding gas. 

From the new AAOmega optical spectroscopy, we have 
found: 

1. Prom 343 spectroscopically confirmed members 
within a 3Mpc radius, the redshift and velocity 
dispersion of Abell 2744 are, respectively, Zdus = 
0.3064 ± 0.0004 and a = 1497 ± 47kms-^ 

2. The velocity distribution is significantly skewed 
and the KMM analysis reveals the velocity distri- 
bution in the central region is characterized by two 
distinct structures in velocity space, separated by 
~ 3000 km s-^ 
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3. Examination of the spatial distribution of the 
galaxies within the two partitions found in the 
KMM analysis of the velocity distribution reveal 
that for the negative peculiar velocity component 
the more significant galaxy concentration lies fur- 
ther to the north, approximately coincident with 
the cool, high mctallicity northern core revealed 
by the Chandra data. 

4. Applying the KMM algorithm to the combined ve- 
locity and spatial information, wc find the statisti- 
cally favored fit is one where the sample is parti- 
tioned into three. Two of the partitions correspond 
to the remnant cores of the two merging clusters in 
the central regions, while the third likely contains 
the galaxies which were in the outskirts of the two 
subclusters and are currently in the process of mix- 
ing to form the new cluster. The third substruc- 
ture also contains those galaxies which belong to 
the northwestern interloper. 

5. The northwestern interloper detected in the Chan- 
dra image is not detected as a dynamical substruc- 
ture, nor as a substructure in local galaxy sur- 
face density when considering only spcctroscopi- 
cally confirmed members. An increase in local 
galaxy surface density is detected when consider- 
ing cluster members defined using photometric red- 
shifts, although only when including galaxies down 
to magnitude R=22.3. 

In combination, the observations presented here have 
been used to outline a skeleton hypothesis for the merger 
history of Abcll 2744. In this scenario, wc have iden- 
tified two significant subclusters in the central region 
which have undergone a violent core passage and are 
now moving away from each other along a roughly north- 
south axis, with a large line-of-sight component. For the 
northwestern subcluster, we propose a scenario where the 
structure is traveling to the north/northeast after peri- 
center. Confirmation of this merger hypothesis requires 
detailed simulations and deeper optical spectroscopy in 
the central regions in order to probe the dynamics of the 
northwestern subcluster. 

The evidence for a post-core passage merger phase re- 
ported here, along with the significant enhancement of 
blue starburst and poststarburst galaxies residing within 
Abell 2744 (Couch et al. 1998), provides further evi- 
dence that major cluster mergers can act as catalysts 
for rapid, environmentally driven evolution in the star 
forming properties of cluster galaxies and that this evo- 
lution occiirs when the chistcr galaxy environment is vi- 
olently rearranged during the core passage phase of a 
major merger. Our large spectroscopic sample, in com- 
bination with deep radio data, can now be used to define 
a large sample of actively transforming galaxies which 
will further test this hypothesis of merger-induced galaxy 
transformation. 
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